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Abstract We review the statistical properties of the main populations of radio sources, 
as emerging from radio and millimeter sky surveys. Recent determinations of local lu- 
minosity functions are presented and compared with earlier estimates still in widespread 
use. A number of unresolved issues are discussed. These include: the (possibly luminosity- 
dependent) decline of source space densities at high redshifts; the possible dichotomies 
between evolutionary properties of low- versus high-luminosity and of flat- versus steep- 
spectrum AGN-powered radio sources; and the nature of sources accounting for the up- 
turn of source counts at sub-mjy levels. It is shown that straightforward extrapolations 
of evolutionary models, accounting for both the far-IR counts and redshift distributions 
of star-forming galaxies, match the radio source counts at flux-density levels of tens 
of pJy remarkably well. We consider the statistical properties of rare but physically 
very interesting classes of sources, such as GHz Peak Spectrum and ADAF/ADIOS 
sources, and radio afterglows of 7-ray bursts. We also discuss the exploitation of large- 
area radio surveys to investigate large scale structure through studies of clustering and 
the Integrated Sachs- Wolfe effect. Finally we briefly describe the potential of the new 
and forthcoming generations of radio telescopes. A compendium of source counts at 
different frequencies is given in an appendix. 
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1 Introduction 



For several decades, extragalactic radio surveys remained the most powerful tool to 
probe the distant universe. Even 'shallow' radio surveys, those of limited radio sensi- 
tivity, reach sources with redshifts predominantly above 0.5. Since the 1960s, the most 
effective method for finding high-z galaxies has been the optical identification of radio 
sources, a situation persisting until the mid-1990s, when the arrival of the new genera- 
tion of 8-10 m class optical/infrared tele scopes, the refurbis hment of the Hubble Space 
Telesco pe, the Lyman-br eak technique i Steidel et alJll996l ) and the Sloan Digital Sky 



survey (|York et alJl200d ) produced an e xplosion of data o n high-redshift galaxies. 



This is not a historical account (see ISullivan 1111120091 ); but listing the revolutions 
in astrophysics and cosmology wrought by radio surveys serves to set out concepts and 
terminology. On the astrophysics side we note the following: 

1. Active Galactic Nu clei (AGNs) The discovery of radio galaxi es feolton et alJ 
ll949l ; lRvle et al.lll950l ') whose apparently prodigious energy release (|Burbidgd l 1959 ) 
suggested Compton catastrophe, calling the cosmological interpretation of redshifts 
into question. 

2. Synchrotron e mission The identification of synchrotron emission (|Ginzburjl95ll : 



IShklovskiill95i ) as the dominant continuum process producing the apparent power- 
law spectra of radio sources. 

Quasars The discovery of quasars, starting with 3C 273 IIHazard et al Jl963l:ISchmidd 



1963), leading to the picture of the collapsed supermassive nucleus i Hovle fc Fowlen 



1963), and hence to the now-accepted view of the pow erful Active Galact ic Nucleus 



(AGN) - mas sive black-hole - accretion disk systems ( Lvnden-B cll 1969) powering 



double-lobed Jjennison fc Das Gupta) | l953f) radio sou rces via 'twin-exhaust' rela- 
tivists beams (|Blandford fc Reeslll974l : IScheuedll974| ). 

Relativisti c beaming The d iscovery of superluminal motions of quasar radio 



comp onents ( Cohen et al.lll97lh . this non-anisot ropic emission (anticipated bv lReed 
1967T ) resolving the Compton non-catastrophe l|WoltieilE966l ) and leading to the 



development of unified models of radio sources: quasars and radio galaxies are one 
and the same, with orientation of the axis to the viewer's line of sight determining 



and the same, with orientation ot the axis to the viewer s line ot sight de termining 
classification via observ ational appearance (|Antonucci fc Miller|[l985l : lBarthellll989l : 
lUrrv fc Padovanj[l995h . 

On the cosmology side we note the following: 

1. Scale of the observable Universe An irrefutable argument bv lRvle fc Scheueil 
( 1955h placed the bulk of 'radio stars' beyond 50 Mpc, and it w as quickly realized 
when arcmin positional accuracy became available (|Smithlll952h that the majority 



of the host galax ies were beyond the reach of the optical telescopes of the epoch. 
iMinkowskil jl960h measured a redshift of 0.46 for 3C 295, the redshift record for a 



galaxy for 10 years. Astronomers had discovered a set of objects substantially 'be- 
yond' thfij^cogjiized Universe. By 1965 the redshift record was 2.0 for the quasar 
3C 9 (Schmidt 1965). Only after the turn of the century did the redshift record be- 
come ^ojrMnelvset by objects discovered in surveys other than at radio wavelengths 
fe.g. lSterijlioOol ). 

History of the Universe Early radio surveys generated a passionate and per- 
sonal debate, the Steady-State vs Big-Bang controversy. It was rooted in the 
simplest statistics to be derived from any survey: the integral source counts, the 
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number of o bjects per unit s k y are a above given intensities or flux densities. As 
discussed by Rvle fc Scheuerl 1 1955h . the source count from the 2C radio survey 



i Shakeshaft et al.lll955 ) showed a cumulative (integral) slope of ~ —3, far steeper 



than that expected from the Steady-State prediction, any reasonable Friedman 
model, or from a static Euclidean universe. For each of these, the initial slope at 
the highest flux densities is —3/2. (Euclidean case: the number of sources, N, is 
proportional to th e volume, i.e. to r 3 f or a sphere; th e flux density is oc r~ 2 , so 
that N oc S~ 3/2 .) iBondi fc Goldl (I1948T ) together with iHovld (|l948h were uncom- 



promising proponents of the new Steady-State theory. Ryle et al. interpreted the 
2C apparent excess of faint sources in terms of the radio sources having far greater 
space density at earlier epochs of the Universe. Confusion, the blending of weak 
sources to produce a continuum of strong sources, was then shown to have disas- 
trous effects o n the early Camb ridge source counts. From an independent survey 
in the South, iMills et al.l ( 19581 ) found an initial slope of -1.65 a fter corr e ctions 



for instrumental effects, significantly lower than that found for 2C. IScheuerl (|l957T l 
developed the P(D) technique, circumventing confusion and showing that the in- 
terferometer results of 2C were consistent with the findings of Mills et al. But the 
damage had been done: cosmologists, led by Hoyle, believed that radio astronomers 
did not know how to interpret their data. 

In 196 5 the 'source-count controversy' became ir relevant in one sen se. IPenzias fc Wilson! 
( 1965) :'ound what was immediately interpreted (|Dicke et al .1 19651 ) as the relic radiation 



from a hot dense phase of the Universe. The Big Bang was confirmed. 

Ryle was right all the time. Integral source-count slopes of —1.8 or even as shallow 
as —1.5 were nowhere near what the known redshifts plus Steady State cosmology - 
or even any standard Friedman cosmology - predicted. These all come out at —1.2 or 
— 1.3, shallower than the asymptotic —1.5 as sources of infinite flux density are not 
observed, and nobody has ever claimed the initial source cou nt slope at any freq uency 
to be as flat as this. The discovery of the fossil radiation fsee IPeebles et al1l2009h may 
indeed have shown that a Big Bang took place; but the source counts demonstrated 
further that objects in the Universe evolve either individually or as a population - a 
concept not fully accepted by the astronomy community until both galaxy sizes and 
star-formation rates were shown to change with epoch. 

Source counts from radio and mm surveys - with errors and biases now under- 
stood - are currently recognized as essential data in delineating the different radio- 
source populations and in defining the cosmology of AGNs. These counts are domi- 
nated down to milli-Jansky (mjy) levels b y the canonical radio s ources, believed to be 
powered by supermassive black-holes (e.g. iBegelman et al.lll984T ) in AGNs. At fainter 
flux-density levels, a flattening of slope in the Euclidean normalized differential counts 
(i.e. co unts of sources with flux density S, within dS, multiplied by S 2 ' 5 . see was 
found ^Windhorst et alJll984 iFomalont et alJll984 ICondon fc Mitchellllliii ). inter- 



preted at the time as the appearance of a new population whose radio emission is, to 
some still-debated extent, associated with star-forming galaxies. 

Radio-source spectra are usually described as power laws (S„ oc the early 

low-frequency meter-wavelength (e.g. 178 MHz) surveys found radio sources with spec- 
tra almost exclusively of steep power-law form, with a ~ 0.8. Later surveys at cm- 



1 We note that this negative sign convention for a is not universal; however the convention 
has been adopted for the K-corrections of optical quasars and for the extrapolation from optical 
to X-rays ('a ox '). 
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Fig. 1 Unified scheme for high radio-power Fanaroff-Riley (1974; FRII) sources (following 
Jackson & Wall 1999). 



wavelengths (higher frequencies, e.g. 5 GHz) found objects of diverse spectral types, 
some with spectra rising to the high frequencies, some with steep low-frequency por- 
tions flattening and rising to the high frequencies, and yet others with a hump in the 
radio regime, or indeed two or more humps. In general, anything which was not 'steep- 
spectrum' in form was called 'flat-spectrum', an inaccurate nomenclature: very few 
truly flat-spectrum sources have been found and even then the flatness persists over 
only a limited frequency range. Nevertheless AGN-powered radio sources are tradition- 
ally classified in two main categories: steep- (a > 0.5) and flat-spectrum (a < 0.5). 
Broadly speaking, to radio telescopes the steep-spectrum objects showed extended 
double-lobed structures, while the flat-spectrum objects were point sources, unresolved 
until the Very-Long-Baseline Interferometry (VLBI) technique provided sub-arc-second 
mapping. The compact nature of flat-spectrum sources led to the conventional inter- 
pretation of synchrotron self-absorption at frequencies below the bump(s), implying 
brightness temperatures of ~ 10 11 K for the estimated magnetic field strengths. 

From a physical point of view, it is appropriate to consider the integrated spectra as 
composites, built of the combination of different components of radio sources. Unified 
models provide a framework for such a discussion. 

In the widely acc e pted 'unificati on' scheme (|Scheuer fe Readheadll979l ; IOrr fe Browne! 
Il982l : IScheuerj|l987l : iBarthell [Tiij ) the appearance of sources, including this steep- 
spectrum/flat-spectrum dichotomy, depends primarily on their their axis orientation 
relative to the observer. This paradig m stems from the discovery of relativistic jets 
( Cohen et al]|l97ll : iMoffet et alJll97a ) giving rise to strongly anisotropic emission. In 
the radio regime (Fig. [TJ , a line-of-sight close to the source jet-axis offers a view of the 
compact, Dopple r-boosted, flat-spectrum base of the approaching jet. Doppler-boosted 
low-radio-power |Fanaroff fc Rilevl Il974l . type I (FRI; edge-dimmed)) sources are as- 
sociated with BL Lac objects, characterized by optically-featureless continua, while 
the powerful type II (FRII; edge-brightened) sources are seen as flat-spectrum radio 
quasars (FSRQs) . The view down the axis offers unobstructed sight of the black-hole - 
accretion disk nucleus at wavelengths from soft X-rays to UV to IR, and this accretion- 
disk radiation may outshine the starlight of the galaxy by 5 magnitudes. The source 
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Fig. 2 Spectral behaviour in the millimeter band of the radio galaxy NGC6251 (left panel) 
and (right panel) 11-GHz isophotes overlaid on the 0.3 GHz map (Mack et al. 1997). The 
low-frequency spectrum is due to the steep-spectrum outer lobes while at higher frequencies 
the flatter-spectrum core-jet system dominates. 



appears stellar, either as a FSRQ or a BL Lac object. FSRQs and BL Lacs are collec- 
tively referred to as blazars. In the case of a side-on view, the observed low-frequency 
emission is dominated by the extended, optically-thin, steep-spectrum components, 
the radio lobes : and the optical counterp art generally appears as an elliptical galaxy. 
A dusty torus 1 Antonucci fe Millerl ll985) hides the black-hole - accretion-disk system 
from our sight (Fig. [TJ . At intermediate angles between the line-of-sight and the jet 
axis, angles at which we can see into the torus but the alignment is not good enough 
to see the Doppler-boosted jet bases, the object appears as a 'steep-spectrum quasar'. 

In general, then, each source has both a compact, flat-spectrum core and extended 
steep-spectrum lobes (Fig. [2}. This already implies that a simple power-law represen- 
tation of the integrated radio spe ctrum can only apply to a limited frequency range. 
The reality is even more complex (Wall 1994). External absorption or, more frequently, 
self-absorption (synchrotron and free-free) can produce spectra rising with frequency 
at the low-frequency optically-thick regime, while at high frequencies the synchrotron 
emission becomes optically thin, pow er law, and energy losses of relativistic electrons 
( "electron ageing" , iKellermannll 196a ) translate into a spectral steepening. 

Two classes of ultra-steep-spectrum (a > 1.3) sources have been discovered. One 
is associated with galaxy clusters; the objects are of relatively low luminosity and gen- 
erally are not associated with any host galaxy. They are diffuse and of several types, 
including cluster 'radio halos', 'radio relics' and 'mini-halos', and each type appears 
to involve reactivation of the hot intra-cluster medium by sh ocks or coolin g flows, the 
observed form depending on the cluster evolutionary state (|Ferettill2o6i ). These 'ra- 
dio ghosts' will not be discussed further here. The second class of ultra-steep-spectrum 
source is very radio-luminous and these are mostly identified with very-high-redshift ra- 
dio galaxies. The high redshifts tempt the suggestion that the steep spectral index is due 
to the effect of redshift moving the steepest part of the spectrum ( where electron ageing 
effects are strong) into the observed frequency range. However, iKlamer et al.l ((2006) 
demonstrated that this is not the dominant mechanism, and that high-redshift radio 
galaxies, discovered by the steep-spectrum technique, have intrinsically power-law spec- 
tra. The selection of ultra-steep-spectrum sources is a very effecti ve, but not the only 
( Jarvis et al.| [2009"). way to find high redshift radio galaxies (see iMilev fc De Breuckl 
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2008, for a comprehensive review), incl uding the one holding t he current record, TN 
J0924-2201 at z = 5.19 (aa| 65 ~ 1.6: Ivan Breueel et al.l l l999l ) . The highest-redshif t 
radio-loud quasar known to date, the z = 6.12 QSO J1427+3312 (jMcGreer et alj|2006l ) . 



also has a steep radio spectrum {a\'\ = 1.1) although it was not discovered through 
this characteristic. 

In very compact regions, synchrotron self-absorption can occur up to very high 
radio frequencies, giving rise to sources with spectral peaks in the GHz range. At high 
radio lumino sities this category comprises the GHz Peaked Spectrum (G PS) sources 
dO'Dealll998h some of which peak at tens of GHz (High Frequency Peakers; lEdge et al.l 
ll998l : lDallacasa et al.ll2000l . l2002l ; lTinti et all2005l ). At low luminosities, high-frequency 
spectral peaks, again due to strong synchrotron self-absorption, may be indicative 
of radiatively inefficient accretion, thought to correspond to late phases of the AGN 
evolution, with luminosities below a few percent of the Edding ton limit (advection- 
dominated accretion flows (ADAF, Quataert fc Nar avanl Il999f l or adiabatic inflow- 



outflow scenarios ( ADIOS, iBlandford fc Begelmanlll999l . 
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As the 'flat' spectra are actually the superposition of emitting regions peakin g 
over a broad frequency range (|Kellermann fc Paulinv-Tothlll969l ; ICotton et alJll980h . 
whose emission is strongly amplified and blue-shifted by relativistic beaming effects, 
a power-law description is a particularly bad approximation. The spectral shapes are 
found to be complicated, and generally show single or multiple humps. Many of these 
show flux-density variations, attributed to the birth and expansion of new components 
and shocks forming in relativistic flows in parsec-scale regions. The variations may 
be on times scales from hours to months or even years, and substantial resources 
have been devoted to monitoring t hese variable sources, led by groups a t Michigan 
(USA) and Metsahovi (Finland) (e.g. lAller et all200ilValtonen et aill2008l '). The latter 
reference shows how global (multi-wavelength and multi-telescope) these monitoring 
programmes have become; moreover the quasi-periodicity for the object in question, 
OJ 287, indicates that it is probably a binary black-hole system. With regard to flux 
variations, we also note the 'Intra-Day Variables' (IDVs), blazars whose flux densities 
vary wildly on time scales from minutes to days: these are flat-spectrum objects with 
extremely small components that show inter-stellar scintil lation (ISS) via the turbulent, 
ionized inter-stellar medium (ISM) of our Galaxy (e.g. lLovell et al1l2007T ). Detailed 
discussion of all these variable objects is beyond the scope of this review. 

T he discovery of Compact Steep Spectrum sources (CSS: lKapahill98ll : |Peacock fc Walll 

Il982l : IO'Dealll998h originally appeared to be an exception to the conventional wisdom 
that steep and flat spectra are associated with extended and compact sources respec- 
tively. CSS sources are unresolved or barely resolved by standard interferometric obser- 
vations (arcsec resolution), and the integrated spectra show peaks at < 0.5 GHz, above 
which the spectral indices (on average, a ~ 0.75) are typical of extended radio sources. 
There is compelling evidence that these objects, as well as GPS and associated types 
of object (HFPs and CSOs - Compact Sy mmetric Objects) are young radio galaxies, 
as summarized concisely bv lSneller] 1 20081 ). 



It follows from the above that the conventional two-population approach (flat- 
and steep-spectrum) assuming power-law spectra is particularly defective at high radio 
frequencies, where several different factors (emergence of compact cores of powerful 
extended sources, steepening by electron energy losses, transition from the optically- 
thick to the optically-thin synchrotron regime of very compact emitting regions, etc.) 
combine to produce complex spectra (see Fig. [3}. Nevertheless, for many practical 
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Fig. 3 Examples of radio-source spectra at mm wavelengths: a flat-spectrum source (top 
left panel); a steep-spectrum source (bottom left panel); a source whose spectrum flattens at 
v ~ 10 GHz (top right panel); a High Frequency Peaker (HF P) source (bottom right panel). 
Data from the NE WPS Catalogue llLopez-Caniego et al.ll2007l . asterisks) and from the AT20G 
Survey ( diamonds, iMassardi et alj|20*08afr 



applications the conventional approach remains useful in describing the bulk population 
properties of AGN-powered radio sources. 

The radio emission of star-forming galaxies is mostly optically-thin synchrotron 
from relativistic electrons interacting with the galactic magnetic field , but with sig- 
nificant free-free contributions from the ionized interstellar medium ICondonl [1992 ; 
iBressan et al . 2002; Clemens et al ] |2008h . At mm wavelengths, however, the radio emis- 
sion is swamped by (thermal) dust emission, whose spectrum rises steeply with increas- 
ing frequency. The w ell-known tight correlation between radio and far-IR emis sion of 
star-forming galaxies l|Helou et al. 1985; Gavazzi et al.lll986l ; ICondon et ai]l99ll ) vastly 
increases the body of data relevant to characterize, or at least constrain the evolution- 
ary properties of this population. However, to date few attempts have been made to 
buil d comprehensive mod els encompassing both radio and far-IR/sub-mm data (but 
sec iGruppioni et alj|2003l ). 

In this paper we first review the observed radio to mm- wave source counts (§0, 
the data on the local luminosity function of different radio source populations (§[3}, 
and the source spectral properties ( §HJ) . Next (§(5| we look at evolutionary models for 
the classical radio sources as well as for individual populations, such as GPS sources, 
ADAF/ADIOS sources, and (§[6]) star- forming galaxies and 7-ray afterglows at ra- 
dio wavelengths. We deal briefly with the Radio Background (§[7} and the Sunyaev- 
Zeldovich effect on cluster and galaxy scales(§[8]). Section [9] contains a summary of the 
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Log(S [Jy]) 

Fig. 4 Differential source counts at 150, 325, 408, 610 MHz normalized to cS„ 2 ' 5 , with c = 
1000, 100, 10, 1 respectively. Reference codes are spelt out in the notes to Tables [TJ [3] and [4] 
The lines are fits yielded by an updated evolution model (Massardi et al., in preparation). 

information on large scale structure stemming from large-area radio surveys. Finally, 
in §[TU]we summarize perspectives for the future, and §QT] contains some conclusions. 

Unless otherwise noted, we adopt a flat ylCDM cosmology with £2^ = 0.7 and 
H = TOkms^Mpc -1 . 



2 Observed source counts 

2.1 From surveys to counts 

The counts of sources demand cosmic evolution, but by themselves provide limited 
information on this evolution. Since even bright radio sources are frequently optically 
faint to invisible, the traditional way to characterize the evolutionary properties re- 
lies heavily on source counts from blind surveys, with limited and incomplete cross- 
waveband identification and redshift information. However, as discovered in the 2C 
survey, getting from a sky survey to a source count is difficult, and modern instru- 
mentation, while generally avoiding the confusion issue which bedevilled 2C, does not 
remove the difficulties. 

It is surface brightness, or rather differential surface brightness above a background 
(CMB, Galactic radiation, ground radiation), which is measured in radio/mm surveys. 
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Fig. 5 Normalized differential source counts at 1.4 GHz. Note that the filled diamonds show 
the counts of AGNs only, while all the other symbols refer to total counts. Reference codes 
are spelt out in the note to Table |51 A straightforward extrapolation of evolutionary models 
fitting the far-IR to mm counts of populations of star-forming (normal late type (spirals or 
sp), starburst (sb), and proto-spheroidal) galaxies, exploiting the well established far-IR/radio 
correlation, naturally accounts for the observed counts below ~ 30 fi Jy (see §[6j . At higher 
flux densities the counts are dominated by radio-loud AGNs: the thick solid line shows the fit 
of the same model as in Fi g. HI The dot-dashed line shows the counts of 7-ray burst (GRB) 
afterglows predicted by the Ciardi & Locb (2000) model. 



Discrete sources stand out from this background by virtue of apparent high differ- 
ential surface brightness, AT^. The simple relations linking ATf, to point-source flux 
density (via the Rayleigh- Jeans approximation and the radiometer equation incorpo- 
rating telescope and rec e iver p arameters) appear in basic radio astronomy texts, e.g 
iBurke fc Graham-Smith! (|l9971 ). 

Surveys are complete only to a given limit in AT^, translating to Jy per beam 
arefl For point sources, this limit is clearly defined. For extended sources, the total 
flux density 

S v = / B v {e,4>)dQ, (1) 

i.e. the incremental brightness B must be integrated over the extent of the source to find 
the total flux density. If a source is extended and its brightness temperature is constant 
across the beam response, then given the Rayleigh- Jeans approximation B = 2krjT) J /}? 

2 1 Jy (Jansky) = 10 -26 W Hz -1 m -2 or 10 -23 ergem -2 s _1 Hz -1 
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Fig. 6 Differential source counts at 4.8 and 8.4 GHz normalized to c X S,7 2 ' 5 , with c = 1 and 
0.1, respectively. Reference codes are spelt out in the notes to Tables |61 and 171 In the upper 
inset sp, sb, and sph stand for spiral, starburst and proto-spheroidal galaxies, respectively. 
The fit to the 4.8 GHz counts is from the same mod el as in Figs. [4] and [5] while at 8.4 GHz 
we show the fit yielded by the Dc Zotti ct al. (2005() model, tailored for dat a above 5 GHz. 
The d otted line shows the counts of GRB afterglows at 3 GHz predicted by the Ciardi & Locb 
(2000) model. 



(kg is the Boltzmann constant, A is wavelength), for a survey sensitivity limit of 5i; m 
per beam, we have from eq. ([TJ 

rp _ ^ 'S'lim /~\ 

b ' min ~ 2k B Jdn ( ) 

The integral is the beam solid angle; for a circular Gaussian beam with full width at 
half maximum of FWHM arcsec, this may be approximated as 2.66 ■ 10 -11 FWHM 2 
sterad. Two iconic sky surveys at 1.4 GHz with th e NRAO Very Larg e Array illustrate 
the brightness limit issue. For the FIRST survey (jBecker et all 1 1993 ) with FWHM = 



5 arcsec and S m j n = 1 mjy, eq. (J2]| gives r m ; n ~ 24 K, while for the NVSS survey 
( Condon et al.|[l99ct ) with a 45 arcsec beam and S lri i n = 3 mjy, T m i n ~ 0.9 K. There 



are significant selection effects which arise as a consequence, most notably the lack of 
sensitivity in FIRST to the majority of spiral galaxies, near and far, as well as to low 
surface brightness features such as ghost or relic radiation. The redeeming features of 
its higher resolution are emphasized below. 
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It is a major undertaking to proceed from a list of deflections in Jy per beam, either 
apparently unresolved, or resolved as regions of emission, to a complete catalogue of 
radio/mm sources. In the first place, there is the surface brightness limitation described 
above; in the compromises of survey design it is critical to decide just what popula- 
tion^) of sources will be incompletely represented. There is the issue of overlap: for 
instance Centaurus A, NGC 5128, the nearest canonical radio galaxy, extends over 9 
degrees of the southern sky; there are many discrete distant sources catalogued within 
the area covered by Cen A. There is also the double nature of radio-galaxy emission: 
this requires that components found as individual detections be 'matched up' or as- 
sembled to find the true flux density of single sources, cores as well as double lobes. 
Moreover many sources show extended regions of lower surface brightness which are 
poorly aligned. If source scale is large enough, pencil-beam or filled aperture telescopes 
are better at finding and mapping these than are aperture-synthesis interferometers. 
The issue of 'missing flux' is notorious for interferometers, because of their limited 
response to the longer wavelengths of the spatial Fourier transform of the brightness 
distribution. 

The difficulties have been brought to sharp focus by the excellent de cision to carry 



out t he two major VLA surveys, FIRST |Becker et alJll995h and NVSS (|Condon et all 
1998), both at 1.4 GHz but differing in resolution by a factor of 9. From these highly 
complementary surveys, t he reality of how diffe rent resolutions affect raw source lists 
may be seen immediately (jfBlake fc Wa ll 2002b]). FIRST and NVSS are far more than 



the sum of the parts. The low resolution of NVSS gains the spiral galaxies and much 
other low-surface-brightness detail not seen in FIRST. The relatively high resolution 
of FIRST can be used to sort out the blends and overlaps in NVSS, and it enables 
direct cross-waveband identifications, a shortcoming of the lower NVSS resolution. 
Used together they can provide samples complete on many criteria; but significant 
effort in examining many individual emission features is still required. 

With regard to surveys using interferometers, the noise level in an interferometric 
image is given by: 



' image 



V^kgTsys I J 
Ar] e r] q 



(3) 



where T sys is the system temperature, A is the antenna surface area, rj e is the aper- 
ture efficiency, the ratio of effective collecting area to surface area, r\ q is the sampling 
efficiency depending on digitization levels and sampling rate, t is the integration time, 
-^basc = N(N — l)/2 is the number of baselines, N is the number of antennas, and Av 
is the bandwidth. [r\ e is generally 0.3 to 0.8, and r\ q 0.7 to 0.9.) The integration time 
per pointing needed to reach a detection limit of say S , i; m = 5<Ti ma g can be straightfor- 
wardly obtained from eq. Q ■ The number of pointings necessary to cover a sky solid 



angle Q s with a telescope field of view FOV iqj 

n p = n s /FOV. (4) 

If the integral counts of sources scale as S - , the number of sources detected in a given 
area scales as t^l 2 . For a given flux density, the number of detections is proportional 



3 This assumes uniform response over the field-of-view. The inevitable non-uniformity across 
the FOV implies an additional factor of ~ 2 for uniform sky coverage. The data from separate 
pointings are combined by squaring the relative response to weight the data by the square of 
the signal-to- noise ratio (SNR). If t he beam is app r oxima ted by a Gaussian, then this process 
effectively halves the beam size; see ICondon et ail Jl998f) . 
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to the surveyed area, i.e. to t. Thus, to maximize the number of detections in a given 
observing time it is necessary to go deeper if /3 > 2 and to survey a larger area if 
ft < 2. The 'narrow and deep' vs. 'wide and shallow' argument for maximizing source 
yield always resolves, at radio frequencies, in favour of the latter, because /3 > 2, 
implying a differential count slope of less than —3, has never been observed at any 
flux-density level. On the ot her hand, very steep counts are observed at m illimeter and 



nux-density level. Un the ot her hand, ver y steep counts are observe 
sub-millimeter wavelengths 1 Austermann et aUl2009l ; ICoppin et al 



Compilation of complete and reliable catalogues, complete samples, almost invari- 
ably involves data at other frequencies. Source-component assembly for example is an 
iterative process which may require cross-waveband identification of the host object, 
galaxy, quasar, etc. The identification process leads on to the construction of complete 
samples, complete at both the survey frequency and at some other wavelength, i.e. in 
optical/IR identifications. Such samples are rare and re quire great observ ational effort. 
One of the best known of these, the '3CRR' sampl e ( Laing et al.lll983T ) is a revised 



version of t he rev ised 3C catalogue ( Bennett] fl962D from the original 3C survey of 



lEdge et all (|l959l ). (The sample is also the most extreme sample of high-power radio 



AGN, and its contents are far from typical of the radio-mm survey popul ation.) The 
proce ss will become easier with large-area optical surveys such as SDSS I York et al.l 
l200(t ) and with the advent of synoptic telescopes such as LSST. 

Given complete samples, then, we can compile source counts. (It should be noted 
that these are frequently constructed by approximations from raw deflection lists, to 
circumvent the labour discussed above. Caveat emptor^) Today the task of checking for 
systematic effects from approximations or statistical procedures is made easier because 
the counts from different survey samples - except for the very deepest ones - overlap 
at various flux-density levels. The counts are usually presented in 'relative differential' 
form, the differential counts dN/dS giving the number of sources per unit area with 
flux density S within dS, subsequently and conveniently normalized to the 'Euclidean' 
form, i.e. multiplied by cS 2 ' 5 , with c being a suitably chosen constant. (A uniform 
source distribution in a static Euclidean universe yields dN/dS oc S~ 2 ' 5 as described 
earlier). A summary of the available source counts at different frequencies is given in 
Tables CKEO] (see also Figs. [5HZj> ■ 

In the case of surveys covering small areas, the field-to- field variations arising from 
the source clustering (sampling variance) furt her adds to the uncertainties. The frac- 
tional variance of the counts is l|PeeblesHl98rj| ): 



with 



where 6 is the angle between the solid angle elements dfi\ and dJ?2 5 and the integrals 
are over the solid angle covered by the survey. 



4 The buyer should also beware of confusing as complete samples (a) lists of sources in 
which large volumes of dat a are assembled from diff erent surveys and different completeness 
algorithms (e.g. PKSCat90, Wright & Otrupcek 1990), and (b) spectral samples, in which flux- 
density measurements at different frequencies are assembled to obtain th e integrated spectra 
of samples of sources not nec essarily selected by survey completeness (e.g. iPaulinv-Toth et al.l 
1966; Kcll ermann et al.lll969ft . 
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The angular correlati on function of NVSS an d FIRST sources (see §[9} is consistent 
with a power-law shape ()Blake fc Wallll2002al lbl: lOverzier et al.ll2003h : 

w(6) ~ l(T 3 (0/degr - 8 , (7) 

for angular separations up to at least 4°. Inserting eq. (0 in eq. (|6]) we get 

a 2 = 2.36 x l(T 3 (fi/deg 2 r - 4 . (8) 

The errors given in Tables HHZ1 include this contribution for surveys over areas < 
25deg 2 . 

Differences between source counts for independent fields are in general far larger 
than these errors imply (Condon 2007). There is little doubt that different calibrations, 
beam corrections and resolution corrections are the dominant if not exclusive culprits. 
Further advances in calibration procedures and characterization of the structures of 
faint sources will be required before sampling variance comes to dominate the errors 
in faint counts of radio sources. 



2.2 Low frequency surveys 

Low-frequency surveys have a long and illustrious (but initially chequered) history, as 
we have mentioned. The most extensive ones, both in terms of area (see also §[9j and 
of depth, are those at ~ 1 GHz and at ~ 5 GHz. The NRAO VLA Sky Survey (NVSS; 
ICondon et alJ l998) covers the sky north of 8 — —40° (82% of the celestial sphere) at 1.4 
GHz, down to ~ 2.5 mjy. It has resolution of 45 arcsec FWHM and the raw catalogue 
contains 1.8 x 10 6 entries. It is comp lemented by the Sydney University Molonglo Sky 
Survey (SUMSS; iMauch et alj|2003h at 0.843 GH z. The survey was co mpleted in 2007 
with the Molonglo Galactic Plane Survey (MGPS: lMurphv et alJl2007l ). and now covers 
the whole sky south of declination —30°. 

The VLA 1-4-GHz FIRST survey (for Faint Images of the Radio Sky at Twenty- 
cm; |Becker^er^l][l99^) is the high-resolution (5 arcsec FWHM) counterpart of NVSS, 
and has yielded accurate (< 1 arcsec rms) radio positions of faint compact sources. 
The new catalog, released in July 2008 (format errors corrected in October 2008), 
covers ~ 8444 deg 2 in the North Galactic cap and 611 deg 2 in the south Galactic cap, 
for a total of 9055 deg 2 yielding a list of ~ 816, 000 objects. Northern and Southern 
areas were both chosen to coincide approximately with the area covered by the SDSS. 
The typical flux density detection threshold of point sources is of about 1 mjy/beam, 
decreasing to 0.75 mjy/beam in the southern Galactic cap equatorial stripe. 

Almost full-sky coverage was also achieved at ~ 5 GHz - albeit to a much higher 
flux-density level - by the combination of the Northern Green Bank G B6 survey with 
the S outhern Parkes-MIT-NRAO (PMN) survey. The GB6 catalog J Gregory et all 



Il996l ^l covers the range 0° < S < 75° down to ~ 18mJy/beam, the FWHM major 
and minor dia meters are of 3' .6 and 3 / .4, respectively. The flux-density limit of the 
PMN catalog ^Griffith fc Wright! [l993l ) is typically ~ 30 mjy/beam but varies with 



declination, which spans the range from —87.5° to +10°; the FWHM is of ~ 4'. 2. 
Other large-area, low-frequency surveys: 

- the VLA Low- Frequency Sky Survey (VLSS; ICohen et al"1l2007l) is a 74-MHz con- 
tinuum survey covering the entire sky North of 5 = —30° to a typical point-source 
detection limit of 0.7 Jy; 
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— the Cambridge 6C survey at 151 MHz (jHales et al J 1 19931 . and references therein) 
covers most of the extragalactic sky above 5 = 30°, but gener ally away from th e 
Galactic plane, with 4.2' x 4.2' esc S resolution. The 7C survey |Hales et alj|2007h . 
at the same frequency, covers a similar region of the sky with higher resolution 
(70" x 70"csc(<5)). A somewhat lower-resolution surve y has been carried ou t in the 
low-declination strip 9h < RA < 16ft, 20° < 5 < 35° (jwaldram et al.lll996h . 

— The 8C survey (|Reeslll990l ; lHales et al.lll995l ) covers the polar cap above 8 — 60° 
at 38 MHz with a typical limiting flux density of about 1 Jy/beam. 

— The W esterbork Northern Sky Survey (WENSS; lRengelink etal]|l997l ; lde Bruvn et"afl 
l2000h covers the 3.14 sr north of S = +30° at 326 MHz with 54" x 54" csc(<5) reso- 
lution in total intensity and linear polarization, to a flux-density limit of approxi- 
mately 18 m Jy/beam. 

For more complete references to low-frequency radio surveys, see Tables [TH?] 



2.3 Deep surveys and sub-mjy counts 



The deepest surveys cover small areas of sky on the scales of the primary beams of 
synthesis telescopes; they are carried out with such telescopes in single long exposures, 
or in nested overlapping sets of such exposures. Because source counts are steep, only 
small survey areas are required to obtain large enough samples of faint sources to be 
statistically significant. 



at , 


5, lmJv (jMitchell & Condon! 1 19851; IWindhorst et al. 


Il985l; iHookins et al.1 Il998l; 


Richards 2000:; Bondi et al.||2003; Cilieei et al. 2003; Hopkins et al.|2003; 


Sevmour et al.l 


2004 


; Huvnh et al. 20051; Prandoni et al.l 120061; iFomalont et all 120061; 


Simpson et al. 


2006 


:lBondi et alj 2007t Ivison et al.l 20071; iBondi et al. 2008; Owen & Morrison 2008). 



The point of inflection was originally interpr et ed as signalling t h e eme rgence of a new 
source population (e.g. ICondonlll984al . Il989l) . IWindhorst et al.l (|l985l ) suggested that 
the majority of sub-mjy radio sour ces are faint b l ue ga laxies, presumably undergoing 
significant star formation (SF), and lDanese et all ( 19871 ) successfully modeled the sub- 



mjy excess counts with evolving starburst galaxies, a model that also described the 
IRAS 60 /im counts. 

More recent data and analyses have confirmed that starburst galaxies are indeed 
a major compon ent of the sub-mjy 1.4 GHz source counts, perhaps dominating be- 
low 3-0.1 mJv feenn et al.lll993l; iRowan-Robinson et al.l 



1993: 



Hopkins et al.lll998l. 

2000| ; ISevmour et al.ll2004 bOQgJ ; IMuxIow et al.ll2005l; IMoss et al ll2007i ; IPadovani et al.1 



200S). However, spectroscopic results by Gruppioni et al.l ( 1999b! ') suggested that early- 



type galaxies were the dominant population at sub-mjy levels. Further, it was recently 
suggested (and modeled) that the flattening of the source counts may be caused by 
radi o-quiet AGN (radio-q uiet quasars and type 2 AGN), rather than star forming galax- 
ies l|Simpson et al.1 2006). Distinct counts for high and low-luminosity radio galaxies 
show that low- luminosity FRI-type galaxies prob ably make a substan tial contribu- 
tion to the counts at 1 mjy and somewhat below ( Gendre fc W all 200°]). Based on a 
co mbination of opti c al an d radio morphology as an identifier for AGN and SF galax- 
ies, IFomalont et al. J2006h suggested th at at most 40% of the sub-mjy radio sources 



are AGNs, while 



Padovani et al.1 l|2007bD indicated that this fraction may be 60-80%. 



Huvnh et al.l 1 20081 ) found that the host galaxy colors and radio-to-optical ratios i 



dicate that low- luminosity (or "radio-quiet") AGN make up a significant proportion 



15 



of the sub-mjy radio population. ISmolcic et al.l ( 20081 ). using a newly developed rest- 
frame-colour based classification in conjunction with the VLA-COSMOS 1.4 GHz ob- 
servations, concluded that the radio population in the flux-density range of ~ 50 fiJy to 
0.7 mjy is a mixture of 30-40% of star forming galaxies and 50-60% of AGN galaxies, 
with a minor contribution (~ 10%) of QSOs. 

The origin of these discrepancies can be traced to three main reasons (see also §[3|). 
First, the identification fraction of radio sources with optical counterparts, which is 
generally taken to be representative of the full radio population, spans a wide range 
(20% to 90%) in literature depending on the depth of both the available radio and 
optical data. Second, it is important to make a distinction between the presence of 
an AGN in t he optical counterpa rt of a radio source, and its contribution to the ra- 
dio emission ( Seymour et aljfeoosj ). Non-radio AGN indicators like optical/IR colours, 
emission lines, mid-IR SEDs, X-ray emission, etc. are not well correlated with the ra- 
dio emission of the AGNs and therefore are not necessarily vali d diagnostics of radi o 
emission powered by accretion onto a supermassive black hole |Mux1ow et all 120051 ). 
Third, there are uncertainties in specifying survey level: deep surveys normally cover 
but one primary beam area, heavily non-uniform in sensitivity. A survey claimed com- 
plete at some specified flux density in the central region alone is in fact heavily biased 
to sources of 5 to 10 times this flux density; the survey as a result is biased to the 
hig her- flux-density popu lation, namely AGNs. 

ISevmour et al.l 1 2008h used four diagnostics (radio morphology, radio spectral index, 
radio/near-IR and mid-IR/radio flux-density ratios) to single out, in a statistical sense, 
radio emission powered by AGN activity. They were able to calculate the source counts 
separately for AGNs and star-forming galaxies. The latter were found to dominate 
below ~ 0.1 mjy at 1.4 GHz, while AGNs still make up around one quarter of the 
co unts at ~ 50 ^Jy- 

iBondi et al.l ( 20081 ) pointed to evidence of a decline of the 1.4 GHz counts below ~ 
0.1 mjy. It is possible that a new upturn may be seen at ;$ 1 pjy, due to the emergence 

of normal star- forming galaxies ( Windhorst et allll999l ; lHopkms et alj |2000). 

Essentially all surveys and catalogues are carried out and compiled without refer- 
ence to polarization (the NVSS being an important exception): linear polarization is 
generally less than a few percent, and certainly at mjy levels, below the uncertain- 
ties in flux densities due to calibration, noise and confusion. An average of circular 
polarizations is generally used. (Subsequent to surveys, thousands of measurements of 
polarization on individual sources have been carried out at different frequencies, with 
the rotation measures thus derive d used to map the details of the Galactic magnetic 
field - see e.g iBrown et al.l (|2007h ). The DRAO 1.4-GHz survey of the ELAIS Nl field 



held - see e.g. Iljrown et alj QzUUYl ) ). Tne UKAU 1.4-liliz survey ol trie JiljAlb INI neld 
( Taylor et al.ll2007h was carried out expressly to examine polarization statistics. The 
data at the faintest flux densities, 0.5 to 1.0 mjy, show a tren d of increasing pola riza- 
tion fraction with decreasing flux density, previously noted bv lMesa et al.l 1 20021 ) and 
iTucci et al. (2004), at variance with current models of population mix and evolution. 



2.4 High frequency surveys and counts 

High-frequency surveys up into the mm-wavelength regime vitally complement their 
low-frequency counterparts. The early cm-wavelength surveys (Parkes 2.7 GHz, NRAO 
5 GHz) in the late 1960s and 1970s found that flat-spectrum sources - or at least 
sources whose integrated spectra were dominated by components showing synchrotron 
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Fig. 7 Differential source counts normalized to S"^ 2 ' 5 for the 15 G Hz 9C survey 
(Wal dram et al J 120031 l2009f) , for the 20 GHz ATCA Bright So urce Sample llMassardi et all 
I2008al) , and for the WMAP 23 GHz survey l lMassardi et al.l2009l) . The model bv lDe Zotti et ail 
(2005) is also shown for comparison (solid line: 20 GHz, dotted: 15 GHz). 

self-absorption - constitute 50% or more of all sources in high flux-density samples. 
Modelling space density to examine evolution demands determination of the extent 
and nature of this emergent population, most members of which are blazars. 

High-frequency surveys are very time-consuming. For telescopes with diffraction- 
limited fields of view the number of pointings necessary to cover a given area scales 
as v 2 . For a given receiver noise and bandwidth, the time per pointing to reach the 
flux level 5* scales as S~ 2 so that, for a typical optically-thin synchrotron spectrum 
(S oc !/ -0 ' 7 ), the survey time scales as However usable bandwidth is roughly 

proportional to frequency, so that the scaling becomes ~ v 2A \ but a 20 GHz survey 
still takes more than ~ 25 times longer than a 5 GHz survey to cover the same area to 
the same flux-density limit. 

High-frequency surveys have an additional aspect of uncertainty: variability. The 
self-absorbed components are frequently unstable, young and rapidly evolving. Vari- 
ability by itself would not be an issue except for the fact that it leads to serious biases. 
This is primarily because a survey will always select objects in a high state at the 
expense of those in a low state, and the steep source count at high flux densities exac- 
erbates this situation. A second issue concerns the spectra. Sources are predominantly 
detected 'high'; to return after the survey for flux-density measurements at other fre- 
quencies guarantees (statistically) that these new measurements will relate to a lower 
state. Non-contemporaneous spectral measurements - if above the survey frequency - 
will be biased in the sense of yielding spectra apparently steeper than at the survey 
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epoch. The bias can have serious consequences for e.g. K-corrections in space-density 
studies, as described below. 

Cosmic Microwave Background (CMB) studies, boosted by the on-going NASA 
WMAP mission and by the forthcoming ESA Planck mission, require an accurate 
characterization of the high-frequency properties of foreground radio sources both in 
total intensity and in polarization. Radio sources are the dominant contaminant of 
small-scale CMB anisotropics at mm wavelengths. This can be seen by recalling that 
the mean contribution of unresolved sources with flux Si to the antenna temperature 
Ta measured within a solid angle Q is: 

Ta = «^tt = L 2 ■ ( 9 ) 

where &b is the Boltzmann constant, A is the observing wavelength, and we have 
taken into account that for high multipoles (£ S> 1), Q ~ (2tt/£) 2 . If sources are 
randomly distributed on the sky, the variance of T4 is equal to the mean, and their 
contribution to the power spectrum of temperature fluctuations grows as I 2 while the 
power spectra of the CMB and of Galactic diffuse emissions decline at large ts (small 
angular scales). Therefore, Poisson fluctuations due to extragalactic sources are the 



dominant contaminant of CMB maps on scales < 30', i.e. t > 400 jPe Zotti et all 

ll999l ; lToffolatti et al.lll999h . 

The diversity and complexity of radio-source spectra, particularly for sources de- 
tected at the higher frequencies, make extrapolations from low frequencies, where ex- 
tensive surveys exist, unreliable for the purpose of establishing CMB contamination. 
Removing this unc ertainty was the pr i mary motivation of th e Ryle- Telescope 9C sur- 



veys at 15.2 GHz (jTavlor et alJbOQll ; IWaldram et al]|2003h . These were specifically 



designed for source subtraction from CMB maps produced by the Very Small Array 
(VS A) at 34 G Hz. The surveys have c overed an area of ~ 520 deg 2 to a ~ 25 mjy com- 
pleteness limit. IWaldram et al.l (|2009h reported on a series of deeper regions, amounting 
to an area of 115 deg 2 complete to approximately 10 mjy, and of 29 deg 2 complete to 
approximately 5.5 mjy. The counts over the full range 5.5 mjy - 1 Jy are well described 
by a simple power-law: 

dN ,1 ( S V 2 ' 15 t -1 -1 tm\ 
— c5l(^— j Jy sr . (10) 

A 20- GHz survey of the full Southern sky to a limit of ~ 50 mjy has been carried out by 
exploiting the Australia Telescope Compact Array (ATCA) fast-scanning capabilities 
(15° min -1 in declination along the meridian) and the 8-GHz bandwidth of an analogue 
correlato r. The correlato r was originally developed for the Taiwanese CMB experiment 
AMiBA l|Lo et al.l lioOll ) but has been applied to three of the 6 22 m dishes of the 



ATCA. A pilot survey ijRicci et al.ll2004l : ISadler et all 120061 ) at 18.5 GHz was carried 
out in 2002 and 2003. It detected 173 sources in the declination range —60° to —70° 
down to 100 mjy. The full survey was begun in 2004 and was completed in 2008. 
More than 5800 sources brighter that 45 mjy were detected below declination 5 — 0° . 
An analysis of a complete flux-limited sub-s ample (5*?nr:H^ > 0-50 J y) comprising 320 



extragalactic radio sources was presented bv lMassardi et aL (2008a) 



Shallow (completeness levels ^ 1 Jy) all-sky surveys at 23, 33, 41, 61, and 94 
GHz have been carried out by the Wilkinson Microwave Anisotropy Probe (W MAP). 
Analyses of WMAP 5-year data have yielded from 388 |Wright et al.l l2009h to 516 
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( Massardi et alJ l2009h detections. Of the latter, 457 are identified with previously- 
catalogued extragalactic sources, 27 with Galactic sources; 32 do not have counterparts 
in lower frequency all sky surveys and may therefore be just high peaks of the highly 
non-Gaussian fluctuation field. 

Co unts at ~ 30 GHz have been estimated from DASI data o ver the range 0.1 t o 
10 Jy |Kovac et alj|2002h . from CBI maps in the range 5-50 mjy |Mason et al.ll2003l ). 
an d down to 1 mjy fro m the SZA blind cluster survey (Muchovej et al. 2009, in prep.). 

IClearv et all (|2005l ) used 33-GHz observations of sources detected at 15 GHz to ex- 
trapolate the 9C counts in the range 20 mjy < S33 < 114 mJv. lMason et all ( 20091 ) car- 
ried out Green Bank Telescope (GBT) and Owens Valley Radio Observatory (OVRO) 
31-GHz observations of 3165 NVSS sources; 15% of them were detected. Under the 
assumption that the SaiGHz/Sl^GHz flux ratio distribution is independent of the 
1.4 GHz flux density over the range of interest, they derived the maximum like- 
lihood 1.4 to 31 GHz spectral index distribution, taking into account 31-GHz up- 
per limits, and exploited it to estimate the 31-GHz source counts at mjy levels: 
N(> S) = (16.7 ± 0.4) deg" 2 (SyimJyr - 80±0 - 01 (0.5 mjy < S < lOmJy). The de - 
rived counts were found to be 15% lower than predicted by the |Pe Zotti et al.l (|2005f ) 
model. 

Preliminary indications of a spectral steepening of flat-spect rum sources above 
~ 20 GHz, beyond th e expectations of the blazar sequence model ( Fossati et al.l fl998; 
iGhisellini et a.1.1 ll99Sl ) have been reported. IWaldram et al.1 1 20071 ) used the spectral- 
index distributions over the range 1.4-43 G Hz based on 'simultaneous' multifrequency 
follow-up observations dBolton et al.| [2004) of a sample of extragalactic sources from 
the 9C survey at 15 GHz to make empirical estimates of the source counts at 22, 
30, 43, 70, and 9 GHz by ext r apola ting the power-law representation of the 15-GHz 
counts [eq. (|10[)]. ISadler et alJ ( 20081 ) carried out simultaneous 20- and 95-GHz flux 
densities measurements for a sample of AT20G sources. The inferred spectral-index 
distribution was used to extrapolate the AT 20G counts to 95 GH z. The extrapolated 
counts are lower than those predicted by the lDe Zotti et all 1 20051 ) mo del, and (except 
at the brightest flux densities) also lower than the extrapolation by iHoldawav et al.1 
dl994Tl of the 5-GHz counts. On the other hand, they are within the range of the 
IWaldram et~aH (120071 ) estimates in the limited flux density ran ge where both data sets 
are valid, althou g h the slopes are significantly different. Both IWaldram et al.l ( 2007f ) 
and lSadler et al.1 ( 20081 ) assume that the spectral index distribution is independent of 
flux density. This can only be true for a limited flux density interval, since the mixture 
of steep-, flat-, and inverted-spectrum sources va ries with flux dens ity. In fact, the 
median 20-95 GHz spectrum (a = 0.39) found b vlSadler et all [l200Sh is much flatter 
than that (a = 0.89) measured at 15-43 GHz by IWaldram et al.l ( 20071 ) for a fainter 
sample. 

Of course, extrapolations from low frequencies can hardly deal with the full com- 
plexity of source spectral and variability properties, and may miss sources with anoma- 
lously inverted spectra falling below the threshold of the low-frequency selection. They 
are therefore no substitute for direct blind high-frequency surveys. On the other hand, 
the recent high frequency surveys (9C, AT20G, WMAP) did not produce "surprises", 
such as a population of sources not present in samples se lected at lower frequencies. 
The analysis of WMAP 5-yr data bv lMassardi et al"T ( 20091 ) has shown that the counts 
at bright flux densities are consistent with a constant spectral index up to 61 GHz, 
although at that frequency there is a marginal indication of a spectral steepening. The 
WMAP counts at 94 GHz are highly uncertain because of the limited number of de- 
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Fig. 8 Local luminosity functions at 1,4 GHz of rad i o AG Ns (red dots) and of star-forming 
galaxies (blue dots), as estimated bv iMauch fc Sadler (2003). The lines show extrapolations to 
1.4 GHz of the 60/^m local luminosity functions of "warm" (usuall y interpreted a s starb urst; 
dashed) and "cold" (normal late type; dotted) IRAS galaxies by Takcuchi ct al. (2003); the 
solid lines are the sum of the two contributions. The yellow lines refer to the linear radio/far-IR, 
relationship of eq. 1141 . while the black lines are based on that of eq. 1151 . which deviates from 
linearity at low luminosities. 



tections and the lack of a reliable flux calibration. How ever, taken at f a ce va lue, the 
WMAP 94-GHz counts are below the predictions by the |Pe Zotti et all (120051 ) model 
by ~ 30 %. This indication is confirmed by recen t measurements of the QUaD collab- 
oration ( Friedman fc QUaD Collaboration! l2009h who suggest that the model counts 



should be rescaled by a factor of 0.7 and of 0.6 at 100 and 150 GHz, respectively. 

An indication in the opposite direction, albeit with very poor statistics, comes 
from the MAMB O 1.2-mm (250 GHz) blank-field imaging survey of ~ 0.75 deg 2 by 
IVoss et al. (2006). This survey has uncovered 3 flat-spectrum radio sources brighter 



than 10 mjy, corresponding to an areal density several times higher than expected 
from extrapolations of low-frequency counts without spectral steepening. 

A 43-GHz survey of ~ 0.5 deg 2 , carried out with ~ 1600 snapshot observations with 
the VLA in D-configuration, found only one certain source down to 10 mjy (Wall et 
al., in preparation). A statistical analy sis of the survey data y ielde d a source-count law 
in good agreement with predictions of fWaldra m et al.l (|2007l ) and lSadler et al] (|2008l ) . 
There is no strong indication of a previously unrecognized population intruding at this 
level. 
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Fig. 9 Contributions of steep- (upper panel) and flat/inverted-spectrum (lower panel) sources 
to the local luminosity functions at 1.4 GHz. In both panels the solid line shows Mauch 
& Sadler's (2007) estimate of the total (i.e., flat/inverted- plus steep-spectrum) local radio 
luminosity function of AGNs. In the upper panel, the dashed line shows Mauch & Sadler's 
(2007) estimate of the local luminosity function of starburst galaxies (mostly steep-spectrum), 
that dominate for log L(ergs _1 Hz - 1 ) < 30; the estimates by Toffolatti et al. (1987) and 
Dunlop & Peacock (1990) of the local luminosity functions of steep-spectrum sources include 
both starburst galaxies and AGNs. The Toffolatti's local luminosity function of fiat-spectrum 
sources (lower panel) joins smoothly with the estimated local luminosity function of BL Lac 
objects obtained by Padovani et al. (2007). The cross-hatched area shows the range spanned by 
the estimates by Pierpaoli &: Perna (2004) of the local luminosity functions of ADAF sources. 



3 Local luminosity functions 



The local luminosity function (LLF) describes the local space density of sources as 
a function of luminosity: it constitutes an important boundary condition for evolu- 
tionary models. Its de termination is complicated by several factors discussed e.g. by 
iToffolatti etHl (|l987T ). Ideally we would like to have a large, complete radio-selected 



sample of sources, all with redshift measurements, all at low enough redshifts for evo- 
lutionary effects to be insignificant, yet distant enough for the redshifts to be accurate 
distant estimators. The sources should be distributed over large enough volumes for 
clustering effects to average away. 

In practice, however, the well-known fact that the redshift distribution of complete 
samples of radio sources peaks at z ~ 1 for all flux-density levels down to ~ 10 mjy, im- 
plies that local sources are swamped by the much more numerous distant ones. Singling 
them out by means of complete redshift measurements is therefore impractical, and 
we must confine ourselves to those brighter than some optical magnitude limit, i.e. we 
must deal with both radio and optical selection. If the magnitude limit for redshift mea- 
surements is too shallow, we lose the contribution of optically-faint galaxies. To some 
extent this bias may be corrected using bivariate (optical/radio) luminosity functions 
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( Feltenlll97i : lAuriemma et alJll977h . Alternatively, radio surveys of optically-selected 
samples can be used. 

Spectacular progress has been recently achieved combining large-area spectroscopic 
surveys (Las Campanas SPSS, 2dF and 6dF) with the NVSS and FIRST surveys 
dMachalski fc Godlowskibood : IMagliocchetti et al.|[20o3 ; ISadler et al1l2002l; best et al ' 



120051: iMauch fc Sadler! 120071 ). A more local sample has been used by I Condon et al 
|2002h . 

If full redshift information is available for a flux-limited radio-selected sample con- 
taining too few local sources for a meaningful LLF to be directly derived, the LLF can 
be estim ated from the lu minosity distribution of the sample, for any chosen evolution 
function (Wal l et al.ll 19801 1 . The results are, by definition, model dependent, although 
the evolution function may be tightly constrained by source counts and other data. 

The key to this process is disentangling the star-forming galaxies from the AGNs. 
Radio AGNs dominate above S1.4GHZ — lOmJy; at lower flux densities an increasing 
fraction of nearby galaxies whose radio emission is fuelled by active star formation ap- 
pears. Optical spectra can be used to identify the dominant process responsible for the 
radio emission of each source. Star- forming galaxies have spectra typical of HII regions 
with strong narrow emission lines of Ha and H/9, while AGNs may have a variety of 
spectra, including pure absorption lines (like spectra of giant elliptical galaxies), LINER 
or conventional type 1 or type 2 AGN spectra. Optical AGN spectra, however, do not 
necessarily imply that the radio emission is of nuclear origin. In fact, there is a body of 
evidence that the star-formation and nuclear activities are tightly connected, but the 
radio and optical emissions of AGNs are largely uncorrelated - about 90% of AGNs 
are radio-quiet. An important diagnostic tool to distinguish between galaxies whose 
radio emission is due to star formation and those harbouring a radio-loud AGN is the 
very-well-established, remarkably tight and nearly lin ear correlation be t ween FIR and 
radio continuum em i ssion from star- forming galaxies ( Helou et al. 19851: 



Il98d : ICondon et al.l fl99ll) . A frequently-used criterion (jCondon et al.l 



Gavazzi et al.l 
2002f ) is that 



galaxies with radio to far-IR flux ratio more than thre e times higher than the mean 
for star-forming galaxies are classified as AGN-powered. IMauch fc Sadlerl \200ti ) found 
disagreement between spectroscopic classification and the radio/FIR diagnostic at the 
~ 10% level; a similar reliability was estimated for their classification based on optical 
spectroscopy. Objects with composite AGN + starburst radio emission are probably 
a primary source of classification ambiguity. This suggests that the classification un- 
certainties may contribute significantly to the overall errors on the local luminosity 
function of each population. Nevertheless, rather accurate estimates of the separate 
1.4-GHz local luminosity functions of AGNs and starburst galaxies are now available 
(see Fig. [SJ . 

With the star-formers disentangled from the radio AGNs, a further dichotomy in the 
local luminosity function evaluation is required. Evolutionary models for radio AGNs 
generally split the total radio AGN local luminosity function into the contributions 
of the steep- and flat-spectrum sources. As discussed in §f3] this is a rather crude, 
but frequently unavoidable, simplification. (Obviously, also the radio spectra of star- 
forming galaxies must be known for evolution models, but the problem is simpler 
because in most cases the spectra are "steep", with mean a ~ 0.7 and a relatively 
narrow dispersion.) 

The 1.4-GHz selection emphasizes steep-spectrum sources, but the flat-spectrum 
sources may be important in some luminosity ranges. The estimates o f separ a ted local 
luminosity functions for the two populations go back to IWall et al.l 1 198ll ). [Peacockl 
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( 19851 ). iToffolatti et~ 



thereafter. iRigbv et al 



19871 ). ISubrahmanva fe Harnettl 1 19871 ). with little progress 
(|2008r ) estimated that the density of steep-spectrum sources 
with logL(1.4GHz)/ergs- 1 Hz- 1 > 32 is ~ (3.0 ± 1.2) x 10" 7 Mpc" 3 . 

Moderate to low-luminosity flat- or inverted-spectrum sources are mostly classified 
as BL Lac objects. Very weak, inverted-spectrum radio sources in the centers of other- 
wise quiescent ellipticals may correspond to late phases of AGN evolution (ADAF or 
ADIOS sour ces, see § 15.5,21) . The obser vational information on this latter population is 
very limited. IPierpaoli fc Pernal <|2004h assumed that their space density equals that of 
elliptical galaxies brighter than L*, and adopted a log-normal luminosity function with 
mean log L(2.7 GHz)/erg s" 1 Hz" 1 in the range 27-28, and dispersion a = 0.25. As il- 
lustrated by Fig. [9] the data on the local luminosity function of flat /inverted-spectrum 
sources already constrain the space density of these sources. 

As for BL Lacs, a serious hindrance in the determination of the luminosity func- 
tion is their essentially featureless spectrum, complicating (or defeating) redshift de- 
termination. However, s everal lines of evidenc e suggest that their luminosity function 
evolves weakly if at all ( Padovani et al.ll2007al ). so that the useful volume for comput- 
ing the local luminosit y function extends up to substantial redshifts. The estimate by 
Pado vani et al ] J2007al) compares well with the LLF of flat-spectrum sources obtained 
by lToffolatti*et"al.l (Il987l ). On the contrary, high luminosity flat-spectrum sources are 
very rare locally and evolve strongly, so that a model independent estimate of the local 
luminosity function is essentially impossible. 



4 Source spectra and evolution 

This discussion of radio spectra is far from exhaustive: it sets out to serve two purposes. 
One is related to the K-corrections, the correction for spectral form which must be 
used to derive luminosities at rest-frame frequencies. Getting these corrections right 
is essential in determining space density. The second issue concerns relating source 
counts at different frequencies, and in particular modelling the poorly-determined high- 
frequency counts from the well-defined low-frequency counts. This limited discussion 
thus ignores some aspects of radio spectral measurements which are critical to the 
astrophysics of radio AGNs, such as variability and monitoring, mentioned briefly in 
the Introduction. 

We have noted that spectra of radio sources are frequently represented as simple 
power-laws, S oc v~ a , with the spectral index, a ~ 0.8 for steep-spectrum sources and 
~ for the flat-spectrum ones. However, all radio galaxies deviate from this simple be- 
haviour. Various physical mechanisms contribute to shaping the emission spectrum. At 
low rest-frame frequencies spectra generally show a sharp decline with decreasing fre- 
quency, attributed to synchrotron self-absorp tion; a low energy cut-off to the spectrum 
of relativistic electrons may also have a role (|Leahv et al.lll989l ). The decline is mostly 
observed at rest frequencies of tens of MHz, but the absorption turnover frequency can 
be orders of magnitude higher, as in GPS and ADAF/ ADIOS sources. 

In the optically-thin regime, the spectral index of synchrotron emission, the domi- 
nant radiation mechanism encountered in classical radio astronomy, reflects the index 
of the energy distribution of relativistic electrons. This distribution is steepened at high 
energies by synchrotron losses as the source radiates, and by inverse Compton losses on 
either the synchrotron photons themselves or on photons of the external environment 
( Krolik fc Che3ll99lh . Inverse Compton losses off the cosmic microwave background 
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(CMB) increase dramatically with redshift since the radiation energy density grows as 
(l + zf. As a consequence, a decline with increasing redshift of the frequency at which 
the spectral steepening occurs can be expected. 

While inverse Compton losses are most important to sources with weak magnetic 
fields, powerful sources may possess more intense magnetic fields enhancing the syn- 
chrotron emission. The faster electron energy losses yield a more pronounced steep- 
ening, correlated with radio power (P). Disentangling the effects of radio power and 
redshift is difficult because in flux-limited samples the more powerful sources are pref- 
erentially found at higher redshifts. A further complication arises because a convex 
spectral shape means that redshifting produces an apparent systematic steepening of 
the spectrum between two fixed observed frequencies as redshift increases. Since the 
redshift information is frequently missing, K-corrections cannot be applied and a P—a 
correlation may arise from any combination of these th ree causes. 

This is the situation for the correlation reported bv lLaing fc Peacock! (|l980l ). Em- 



ploying (a large proportion o f) redshift estimates for a sample drawn from the 38 MHz 
8C survey. FLacv et all 1 19931) showed that the high-frequency (2 GHz) spectral index 



correlated more closely with redshift than with luminosity. While at first sight this may 
suggest the domin ant importance of inverse Compton losses on high-frequency spectra, 
lLacv et all (jl993T ) pointed out that the correlation between spectral index and redshift 



weakens when the radio K-correction is applied. This means that such correlation may 
be induced, at least in part, by the spectral curvature due to self absorption at the 
very low selection frequency. In fact, magnetic fields in the very luminous Lacy et al. 
source s should be strong en ough for synchrotron losses to dominate inverse Compton 
losses, blunder! et"aD l| 19991 ). studying a number of complete samples of radio sources 



selected at frequencies close to 151 MHz, with a coverage of the P—z plane (see § 
substantially improved over previous studies, concluded that: 

1. The rest-frame spectral index at low frequency depends on the source luminosity 
(P—a correlation), but also on physical size (D-a correlation) in the sense that 
sources with larger physical sizes D have steeper spectra. 

2. The rest-frame spectral index at high frequency (GHz) depends on the source red- 
shift. 

Simple expressions for the average rest-frame spectra of FRI and FRII radio ga laxies 
as a funct ion of radio pow er and Fanaroff-Riley type (|Fanaroff fc Rilevl 11974 ) were 
derived by I Jackson fc Walll (|200lh . 

With regard to the second issue, relating source counts at different frequencies, 
relevant aspects are to what extent a power-law approximation of the source spectra 
may be viable, i.e. to what extent low-frequency self-absorption, electron ageing effects 
at high frequencies etc. can be neglected; up to what frequencies do blazars have "flat" 
spectra; are source spectra correlated with other parameters (luminosity, redshift) and, 
if so, how can these correlations be described? 

In this regard we have noted that surveys at frequencies of 5 GHz and higher 
are dominated (at least at the higher flux densities) by 'flat-spectrum' sources. The 
spectra of these sources are generally not power-laws, but have complex and individual 
behaviour, showing spectral bumps, flattenings or inversions (i.e. flux increasing with 
increasing frequency), frequently bending to steeper power-laws at higher frequencies. 
Examples are shown in Fig. [3] The complex behaviour results from the superposition 
of the peaked (self-absorption) spectra of up to several components. These components 
are generally beamed relativistically with the object-axis close to the line of sight; they 
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are the parts of jet-base components racing towards the observer at highly relativistic 
speeds. 

The dominant populations of flat-spectrum sources are BL Lac objects and flat- 
spectrum radio quasars (FSRQs), collectively referred to as "blazars". Their spectral 
energy distributions (i.e. the distributions of uL v ) show two broad peaks. The low- 
energy one, extending from the radio to the UV and sometimes also to X-rays, is 
attributed to synchrotron emission from a relativistic jet, while the high-energy one, 
in the 7-ray range, is interpreted as an inverse-Compton component arising from up- 
scattering of eith er the synchrotron photons themselves (sync hrotron self-Compton 
process, SSC, e.g. lMaraschi et al.lll992l ; lBloom fc Marscheril 19931 ) or photons produced 
by the accretion disc near the central black-hole and /or scattered/reprocessed in the 
bro ad-line regi o n felandfordl 1 19931 : ISikora et al] Il994h . 

iGear et al.1 l| 19941 ) investigated the radio to sub-millimeter spectra of a random 
sample of very luminous BL Lacs and radio-loud violently variable quasars. They found 
generally flat or slowly rising 5-37 GHz spectra (median q| 7 ~ for both populations), 
and declining 150-375 GHz spectra, with a statistically significant difference between 
BL Lacs and quasars, the former having flatter spectra (median afrjg — 0.43) than the 
latter (median af 7 ^ ~ 0.73). 

Indication of strong spectral curvature was reported bv ljarvis fc Rawlingsl <|200Ch 
for a qu asi-complete sample draw n from the 2.7 GHz Parkes half-Jansky flat-spectrum 
sample ( Drinkwater et alii 19971 ). It should be noted however that the data used to 
constru ct the radio spectr a are heterogeneous, and the bias which this introduces is very 
serious (|Wall et al.ll2005l ). An extrapolation of this spectral behaviour to frequencies 
above 20 GHz would be in conflict with WMAP finding of a median spectral index 
a ~ dBennett et al.ll2003l ). 

iFossati et"alT "i 19981 ) found evidence for an anticorrelation between the frequency of 
the synchrotron peak v v and the blazar luminosity, and proposed a scenario, dubbed 
"the blazar sequence" , for a unified physical unde rstanding of the vast ra nge of blazar 
properties. The scenario was further extended by iGhisellini et all l| 19981 ), The idea is 
that blazars indeed constitute a spectral sequence, the source power being the fun- 
damental parameter. The more luminous sources are "redder", i.e. have both the 
synchrotron and the inverse Compton components peaking at lower frequencies than 
the lower-luminosity "blue" blazars. If so, the sub-mm steepening could be a prop- 
erty of only t he brightest sources. The validity of the s c enario has been questioned , 
however (e.g., Giommi et al.l 19991; Padovani et al.l 2003]; Caccianiga fc Marchal 12004 ; 



Anton fc Browne82005l;lNieppola et alj2006l;lLandt et alBoOG; Padovani 2007:; lNieppola 
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On the whole, the spectral curvature question is still subject to dispute. Sub- 
stantial progress is expected fro m the forthcoming surveys by the Planck mission 
(|The Planck Collaboration! 120061) . covering the range 30-857 GHz, that will provide 
the first complete samples allowing unbiased studies of the high frequency behaviour. 



5 Evolutionary models: radio AGNs 

By 1966 the counts produced from the Cambridge 3C and 4C surveys at 178 MHz 
( Gower|[l966l ) had already defined the characteristic shape confirmed and refined by 
later surveys: in the Euclidean normalized presentation, a 'Euclidean' portion at the 
brightest flux densities, followed at lower flux densities by a 'bump' (dubbed 'bulge' by 
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IWalj[l994h . and then a roll-off toward the faint intensities. The 'bump' is the signature 
of strong cosmic evolution, implying that the universe must have evolved from a state 
of 'violent activity' in the past to a more quiescent ph ase at the present epo ch, in con- 



tradiction to the steady state model as pointed out bv lRvle fc Clarkg ( 1961 ). It is also 
in co mplete contradiction to uniformly-populated Friedman models (e.g. Wal et al.l 



198d ). A further implication of strong evolution, also noted by iRvle fc Clarke! l|l96ll ) 



is that it swamps the effect on the source counts of different cosmological models, 
frustrating the original hopes that the counts could inform on the geometry of the 
universe. The evolution is indeed spectacular, yielding co-moving source densities at 
z ~ 1-2 a factor of ^ 10 3 higher than those of local sources of similarly high lumi- 
nosities (|Longairill966h . As stressed bv lCondonl(|l989ft , this means that the fraction of 
nearby radio sources is low even at large flux densities; the median redshifts of radio 
sources in complete samples with limits ranging from mjy to Jy are close to 0.8. 

Two reference evolutionary scenarios for radio sources (not to be interpreted liter- 
ally) were used in the 1960s to interpret the counts; they remain popular today. In the 
luminosity evolution scenario, the comoving density of sources is constant but lumi- 
nosities vary with cosmic epoch, while in the density evolution scenario the comoving 
density of sources of any luminosity varies. A density and luminosity evolution of the 
luminosity function at the frequency v can be described as: 

p(L, z; u) = f d (z)p{LISl(*)> * = 0; v), (11) 

where fd(z) and fi(z) are the functions describing the density and the luminosity 

evolution, respectively. 

In his pioneering investigation. lLongai r (1966) showed that either fd(z) = {l + z) m 



or fi(z) — (1 + 2:)™ would fit his data, provided that only sources with high radio lumi- 
nosity (log(Z/]78MHz) ~ 33.5 in the case of luminosity evolution or log(Li78]yiHz) ~ 34.9 

in the case of density evolution; luminosities in ergs -1 Hz -1 ) evolve. This is called 
power-law evolution in which the function is a power of the expansion factor of the 
Universe. (To be precise, Longair quoted evolution functions proportional to t n in an 
Einstein-de Sitter universe, where the cosmic time t is proportional to (l + z)~ 3 / 2 ). The 
power-law evolution diverges at high redshifts, and must be trun cated. The problem 
is avo ided in t he exponential ev o lution model proposed both by iDoroshkevich et al.l 
1 1970ft and by iRowan-Robinsonl (|l970ft . f d (z) = exp(fc-r(z)) or fi(z) = exp(fc'r(z)), 



where r is the lookback time in units of the Hubble time Hq , and k~ 

the evolutionary timescales in the same units. 

A further step in ex ploring radio- AGN evolution was taken bv lSchmidtl (|l968h and 



A tu rtncr step in exploring radio- AOIN evolution was taken by bcnmidt ( 1908 ) and 
Rowan- Robinson! ( 19681 ): the development of the Knax or 'luminosity - volume' test. 



The test avoids mentioning source counts, at the time still lingering under a cloud of 
suspicion. It requires a complete sample, and it requires that if t wo or mo r e limi ts to this 
completeness are in play, they all be considered. For example [Schmidd (l968) applied 



the test to the quasars of the 3CR catalogue, a sample defined by the radio flux-density 
survey limit and by the optical magnitude limit of the original Palomar Observatory Sky 
Survey. The test is simple: take each object and "push" it out in increasing redshift until 
it becomes faint enough to encounter one of the two limits, radio or optical. (The process 
requires that the optical and radio spectra be known so as to define K-corrections in 
each band.) Stop "pushing" at the first limit encountered; this defines z m ax and hence 
Vmax, the maximum co-moving volume in which the object could be found. Calculate 
the co-moving V for the object, the volume defined by the object's redshift, and then 
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form the ratio V/V ma ,x- Do this for all objects in the sample. It is easy to show that 
if they are uniformly distributed in space, the values should be uniformly distributed 
between and 1.0, i.e. (V/Vmax) = 0.5 ± 1/V127V where iV is the number of objects 
in the sample. The statistic is a maximum-likelihood estimator and is unbiased. It 
remains in w i despr ead use and has undergone many refine ments, e.g the C~ me thod of 
iLvnden-Belll (|l97ll ). the application to co mbined samples dAvni fe Bahcanlll980l ). V m 
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methods to extract an evolution function l|Choloniewski|[l987l. and references therein) 



methods to estimate the l uminosity function direct ly (Felten 1976; Eales 19931) . and the 



bande d V/Vmax method ( Dunlop fc Peacocklll99oT ). There is a vast literature; IWillmer] 



(Il997h is a good st arting point. 

ISchmidJ l| 19681 ) foun d strong evolut i on for the 3CR quasars, i n accord with the 
evolution derived by e.g. lLongairl l| 19661 ). lLongair fc Scheueil |l97Cl) showed how very 



closely the test was related to the source-count test and emphasized that the two 
methods were far from independent. However, V/Vmax is versatile, simple in concept, 
model-free, and comes with a statistical pedigree, in contrast to source counts, the trial 
and error process of guessing suitable evolution functions, and the frequently less-than- 
transparent statistical methods used to compare these with source-count and redshift 
data. 

However, the evolution functions mentioned ab ove may carry physica l significance 



For example a physical interpretation was offered bv lCavaliere et all 1 1985h who pointed 
out that models for gravitational energy release near collapsed massive objects yield 
dL/dt — —A(t)L 1+p . If A = constant and p = we get the exponential luminosity 
evolution function with timescale fc" 1 = A" 1 . If A(t) oc t" 1 we get a power-law 
evolution. It must be stressed that these evolutionary laws apply to source populations, 
not to individual sources: the evolutionary timescales, fc -1 , are found to be of order 
of a few to several Gyrs, while the lifetimes of brig ht radio AGN phases are at least 
one an d probably two orders of magnitude shorter (|Bird et alJbOQsl ) . iGrueff fc Vigottl 
( 19771 ) made an early attempt to build a model explicitly constrained by astrophysical 



factors. They linked the radio-source formation to that of parent galaxies and assumed 
radio-emitting lifetimes inversely proportional to radio luminosities. 

Most evolution models have ignored the distributions of spectral indices around the 
mean values (w hich may be lu minosity dependent) for both steep- and flat-spectrum 
sources (but see lCondonl 1984al ) . If such distributions can be approximated as Gaussians 
with dispersion o, and the differential counts have a power-law shape, n(S) = fcS'"' 5 , 
the mean spectra l index a of sources with given value o f S shifts with frequency 
(|Kellermannlll964l ; I Condon! [l984al ; [5 ancsc fc de Zottjfl984l ): 



a = a + o- 2 (l-p)ln(v/v ), (12) 

and the amplitude of the counts scales with frequency as fc = fco(f /vo) q , with: 

q = Q (l - p) + 0.5a 2 (l - /3f In(i//wo). (13) 

The dispersion of the mean spectral indices thus counteracts the effect on counts of the 
high-frequency steepening (§2}. The effect is amplified by the increase with frequ ency 



of the variability amplitude ( Impev fc Neugebauerl 1 19881 : ICiaramella et all [2004) . We 



note however that (a) power-law approximations for source counts hold over very lim- 
ited flux-density ranges only; and (b) a Gaussian approximation for spectral-index 
distributions holds only for low-frequency surveys. By 1.4 GHz, spectral-index distri- 
butions have pronounced tails towards flat spectra and by 5 GHz, the spectral-index 
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distribution is almost the sum of two Gaussians with peaks at ~ 0.8 (steep-spectrum 
sources) and ~ 0.0 (flat-spectrum sources). Clearly, the above equations can be used 
to calculate a and q for each population. 



5.1 Low versus high luminosities 

Extensive discussions can be found in the literature on whether or not the cosmological 
evolution is a property of powerful radio sources only. The origin of the discussion is 
lLongairl |l96d )'s classic study. He adopted a luminosity function extending over about 



8 decades in luminosity and with a shape not far from a power-law, so that density 
and luminosity evolution were essentially equivalent over the flux density range covered 
by the counts available at the time. Under these assumptions, the relative narrowness 
of the evolution bump in normalized counts (width less than 2 orders of magnitude) 
compared to the breadth of the luminosity function (8 orders of magnitude) could 
be accounte d for if o nly th e brightest sources evolve strongly. It was a fundamental 
discovery of lLongairl (|l966h 's investigation that 'differential evolution' was therefore 



essential: if all sources evolved equally, there are far too many faint sources, i.e. the 

model bump is far too broad. 

As mo re data accumulat ed, the evolutionary models were refined. iRobertsonl (| 19781 . 

ll98Ch and I Wall et alj (|l980h factorized the evolving luminosity function as p(L, 



F(L, z)po(L), where F(L, z) is the evolution function, able to represent density or lumi- 
nosity evolution, or a combination of both. Models assuming a luminosity-independent 
evolution function were found to produce unsatisfactory results, while (following the 
qualitative description of the previous paragraph) good fits of th e data were obtained 
only by assuming a far stronger F(L, z) for more luminous sources, banese et alJ (1983) 



proposed a luminosity-dependent luminosity evolution model in which the luminos- 
ity evolution timescale increases with decreasing luminosity and exceeds H^ 1 , so that 
sources evolve weakly (if at all), below some critical luminosity L s ~ 10 31 ergs -1 Hz . 
This is in keeping with the indications that a variety of physical processes can sus- 
tain a ste ady low-level fueling of the central engine for ti mes longer than the Hub- 



ble time dCavaliere etahl 119851 ). The [5 anese et al . (1987) model was exploited by 



iToffolatti et all (|l998T l~to carry out remarkably successful predictions of radio source 
contributions to small scale anisotropics measured by Cosmic Microwave Background 
experiments. 

The differential evolution is suggestive of two populations. Pushing the high/low 
luminosity dichotomy to the extreme, some investigators explicitly considered two pop- 
ulations, one of non-evolvi ng low-lumin osity sources, and the other of high-luminosity, 
strongly evolving sources. IWalll (|l98d ) identified the two populations with FRI and 
FRII radio galaxies. The border between FRI and FRII is, approximately, at -Li.4gHz ~ 
10 32 ergs -1 Hz -1 , although the division appe ars to be dependent on both radio power 
and o ptical luminosity of the host galaxy (cf. iLedlow fc Owenlll996l ). Ijackson fc Walll 
( 19991 ) extended the scheme to flat-spectrum sources, assuming that BL Lac objects 



and flat-spectrum radio quasars are the beamed counterparts of FRI and FRII objects 
respectively, as discussed in § 15.41 

Before the advent of 2dF and SDSS sky surveys, the space distribution o f low - 
luminosity radio AGN was a matter of speculation. For example, lLaing et al.1 |l983) 
showed that the low-luminosity radio galaxies of 3CRR [log^^gMHz/erg s _1 Hz _1 )< 
34] gave values of (V/Vmax) ~ 0.50, suggesting little or no evolution. This view 
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was supported by IClewlev fc Jarvisl <|200d ) who found that the comoving space den- 
sity of sources fainter than Li.4gHz ~ 4 ■ 10 32 ergs -1 Hz -1 remains approximately 
constant with increasing redshift up to z ~ 0.5. However, numbers were small and 
uncertainties undoubtedly permitted some mild evolutio n. On the basis of dis cover- 
ing two distant FRI galaxies in the Hubble Deep Field, ISnellen fc Bestl ( 200 lh pro- 
posed that FRI galaxies showed significant evolution. With recourse to the 2dF sam- 
ple of galaxy redshifts, ISadler et al.l (|2007h quantified this: they found substantial 
cosmological evolution over the redshift range < z < 0.7 of radio galaxies with 
10 31 < -t/i.4GHz < 10 32 ergs -1 Hz -1 , i.e. in the luminosity range of FRI sources. 
They also found indications that the most powerful radio galaxies in their sample 
(with I/i.4GHz > 10 33 ergs -1 Hz -1 ) undergo more rapid ev olution over the sam e red- 
shift range. The latter findings are consistent with those bv lWillott et al.l ( 200 ll ). who 
used emission- line strength rather than morphology (FRI/FRII) to discriminate be- 
tween radio-source populations. The critical radio luminosity dividing sources with 
weak/absent emission lines (the less radio- luminous population of FRI and FRII), and 
the more radio-luminous population of strong-line FR II radio galaxies and quasars, is 
approximately Li.4GHz ~ 6 ■ 10 32 ergs -1 Hz -1 . Both populations show evidence for 
evolution, but the comoving density of the more powerful sources rises far more dra- 
matically than that of the low-luminosity population. 

It must be noted, however, that a luminosity-dependent evolution function, F(L, z), 
does not necessarily imply a luminosity-dependent luminosity evolution. For example, 
in the case of a luminosity function p oc L~ 7 a uniform power-law evolution, L(z) — 
L(0)(1 + z)°, yields p[L(z)] = p[L(0)](l + z) al ; the evolution function depends on 
the slope of the luminosity function. If the luminosity function levels off below some 
'bending' luminosity Lj,, a luminosity-independent luminosity evolution translates into 
a constant comoving space density at low luminosities, with strong variati ons with 
epoch c onfined to the steep portion of the luminosity function. As shown by ICondonl 
(1984b), it is then possible to fit a wide range of radio data with a model assuming 
that all sources evolve equally. 

Finally, we note that the debate on the evolution of low luminosity radio AGNs has 
been somewhat muddled for some time by the poor knowledge of their local luminosity 
function. As discussed in §[3l it is now clear that the faint portion of the radio luminosity 
function is dominated by starburst galaxies, while the luminosity function of radio 
AGNs flattens below logL(1.4GHz)/ergs -1 Hz -1 ~ 31. 



5.2 Steep- versus flat-spectrum sources 

The width of the bump in the Euclidean normalized counts increases with increasing 
frequency. The bump of the steep-spectrum population dominating the low-frequency 
counts shifts to fainter flux densities and gradually fades as survey frequency increases. 
But as survey frequency increases the evolution bump of flat-spectrum sources be- 
comes steadily more prominent, combining with the steep-spectrum bump to produce 
an increasingly broad overall maximu m, a broader range of flux densities over which 
the count slope is close to Euclidean (|Kellermann fc Walll Il987h . This initially mis- 
led investigators to consider different space distributions and evolutionary laws for 
steep- and flat-spectrum sources, with less evolution for the latter. Demonstrating the 
lack of independence of V/Vmax. and source-count results, si milar indi c ation s of little 
evolution for flat-spectrum objects came from this direction. ISchmidtl (1968) found a 
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slight indication in the 3CR sample of quasars that (V/Vmax) appeared to be less for 
the flatter-spectrum objects. When substantial flat-spectrum quasar samples became 
available from cm-wavelengt h surveys, it was found that (V/Vma x) ~ 0.5, consistent 
with a uniform distribution ISchmiddl 19761 ; iMasson fc Wall ll977T ) . This is in contrast 
to (V/Vmax) — 0.7 for steep-spectrum quasars, indicative of a strong increase of their 
spa ce density with z . 

IPeacock fc Gulll ( 198ll ) stressed that the data available at the time poorly con- 
stra ined the luminosi ty function over most areas of the luminosity-redshift plane (see 
also IWall et all Il98lh , and pointed out that within the regions where the luminosity 
function was reasonably well defined, steep- and flat-spectrum sources behave similarly: 
both spectral types undergo differen tial evol u tion, t he strength of evolution increas- 
ing with luminosity. The analysis by ICond'onl (1984a]) confirmed that the two spectral 
classes may indeed evolve similarly within the constraints of the data. This position 
received further support fr om comparative ana lysis of evolution of AGNs in the radio, 
optical and X-ray bands bv lDanese et al.l (|l985h . who noted that the apparently weaker 
evolution of flat-spectrum sources might be related to relativistic beaming effects that 
boost the apparent radio luminosity. In particular, BL Lac objects are probably associ- 
ated with low-luminosity steep-spectrum radio sources which also show weak evolution . 
More extensive data and further analyses ()Danese et al.lfl987l : lDunlop fc PeacockfT990h 
reconciled the epoch-dependent spatial distributions of the two populations, showing 
that they were essentially identical. This finding produced a necessary consistency for 
the success of the unified scheme for radio-loud AGNs; it is clearly not possible for side- 
on and end-on versions of the same populations to have different space distributions. 
'Unified evolutionary schemes', in which flat-spectrum quasars an d BL-Lac objects 
are 'b eame d' versions of FR I and FR II sources, were presented bv lUrrv fc Padovanil 
( 19951 ) and I Jackson fc Wall l\199$ ) and are discussed in g [5T4l 

Some of the apparent discrepancy between these analyses can perhaps be resolved 
noting that, while (V/Vmax) > 0.5 means an increase with redshift of the source density, 
(V/Vmax) — 0.5 does not necessarily mean no evolution. Evolution increasing the 
source density up to z ~ 2 and decreasing it afterwards may yield (V/Vmax) — 0.5 for 
quasars visible up to high z; at high-z the visibility of flat/inverted-spectrum quasars 
is enhanced compared to the steep-spectrum ones by the more favourable K-correction. 
This may be part of the explanation; but probably the majority of the explanation for 
disparate (V/Vmax) results arises from selecting the steep-spectrum samples from their 
initial steep portion of source count, while selecting the flat-spectrum sources from an 
effectively fainter (and flatter) portion of the source count. 



5.3 High-z evolution 

At some epoch, radio galaxies and quasars have to be born; some epoch after recom- 
bination has to have assembled galaxies suitable to harbour both massive black hole 
cores and fuel systems for these to produce collimated twin-beam radio AGNs. We have 
described how easily we can trace the very strong increase of the co-moving number 
density of powerful radio sources between redshift ~ and 2. How hard is it to find 
the more distant epoch at which this relatively high co-moving density falls to signify 
the epoch of AGN birth? 

This is an important issue. The radio activity is associated with processes driving 
the growth of supermassive black holes (SMBHs) in galactic nuclei, which in turn is 
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tied to galaxy formation and evolution. The radio activity may also drive feedback 
processes that have a key role in the evolution of black holes and their host galaxies 
( Granato et "all |2004| ; best et al.l [ioo^ ; bower et al]|2006l ; ICroton et ailbood). Power- 
ful radio galaxies a nd quasars trace the most massive SMBHs (jPunlop et al.l 120031: 
McLure et al.ll2004h which are hosted by the most massive galaxies ( Ferrarese fc Fordl 



2005). Radio survey data are unaffected by dust extinction and may thus provide cru- 



cial tests of the optical/IR results on galaxy for mation. 

One of the achievements of the outstanding founlop fc Peacockl l|l990h free- form 
evolution analysis was to provide evidence for a decline (or high redshift "cut-off") in 
the number density of sources with both steep and flat spectra at redshifts beyond 2.5 
to 3.0. However their samples were incomplete in redshift information, and the results 
were dependent on the accuracy of photometric redshifts ascribed to sources of the 

optically-faintest host galaxies. 

Tracing the redshift cutoff is still under investigation, [shaver et all ( 19961 ) reported 
evidence of a strong decrease of the space dens ity of flat-spectrum radio-loud quasars 
(FSRQs) for z > 3. J arvis fc Rawlingsl l|2000h disputed this, pointing out that the 
apparently curved nature of many of the radio spectra involved, and in particular a 
spectral steepening to the high frequencies, might reduce or remove the significance of 
an apparent redshift cutoff. Ignoring such steepening leads to the prediction of more 
high-z sources than are actually seen, and this could be misinterpreted as evidence for 
a decline of the space density. They concluded that the comoving volume covered by 
the available samples is too small to make definitive statements about any redshift cut- 
off for the most luminous flat-spectrum sources, alth ough both a c onsta nt comoving 
density and a decline as abrupt as those envisaged bv lShaver et~ai1 [| 19961 ) were found 
to be ma rginally inconsiste nt with the data. Using a larger sample with a more rigorous 
analysis. IWall et al.l 1 20051 ) found evidence, significant at > 4<r, of a diminution in space 
density of flat-spectrum q uasars at z > 3, consistent with the redshift-cutoff forms 
of both optically-selected dSchmidt et aL 19951 ; iFan et all |2004| ) and X-ray selected 



( Hasinger et alj 20051 ; ISilverman et al. 20051 ) quasars. Wall et al. drew attention to a 



major source of bias in the Jarvis et al. sample, the non-contemporaneous nature of the 
radio flux-density measurements. It is certain that this introduces too much apparent 
high-frequency curvature and decline into the radio spectra. The effect of the bias was 
demonstrated by Wall et al., using both contemporaneous and non-contemporaneous 
spectral data. 

Using three samples selected at low frequencies. I Jarvis et alj (2001) found that the 
space density of the most radio-luminous steep-spectrum radio sources is consistent 
with being c onstant between z ~ 2.5 a nd z ~ 4.5 and excluded a decline as steep as 
sugges ted by IShaver et al. I l| 19961 .1 19991 ). This conclusion was confirmed bv lCruz et alj 
(2007). However, the samples re main small and incomplete in redshift information and 
as for lDunlop fc Peacockl ( 1990l ). the faintest host galaxies require redshift estimation 
from a K — z plot. It is much harder to track any redshift diminution for steep-spectrum 
radio galaxies then for FSRQ and the statistical uncertainties are inevitably greater; 
showing that the steep-spectrum samples are consistent with a uniform distribution 
does not disprove the redshift cutoff found by I Wall et al.l (2005) for FSRQ. 

These apparently contradictory resu lts may again be accounted for by luminosity- 
dependent evolution. iHook et al. I |l998l ) reported indications that the high- z decline of 
the space density of flat-spectrum quasars is more pronounced and starts at lower red- 
shifts for less powerful sources. These indications were confirmed by subsequent studies 
that did not distinguish between flat- and steep-spectrum sources ( Waddington et all 
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Fig. 10 Relative space density of flat-spectrum radio quasars from the Par kes Quarter-Jansk y 
sample, as derived from a Maximum-Likelihood analysis similar to that of I Wall et al. (2008). 
The s ample and Single-Object Survey (SOS) technique used were described by IWall et alj 
J2005f) . The peak of FSRQ activity is a monotonic function of 2.7-GHz monochromatic radio 
luminosity (L in units of ergss -1 Hz -1 ) in the pseudo-downsizing sense, i.e. lower luminosities 
have the peak of their activity at lower redshifts. The green curve in each diagram represents the 
global solution for the entire sample. This 'down-sizing' is similar to that found for quasars 
selected at X-ray wavelengths (Hasinger et al. 2005) and for submm galaxies (SMG) from 
JCMT surveys (Wall, Pope & Scott 2008). 



l200ll ; IVigotti et al.ll2003l ; ICirasuolo et al.ll2005l . I2006I ). The luminosity dependence of 
the high- z decline is qualitatively s imilar to the do wnsizinq observed for galaxies and op- 
ticall y and X-ray selected quasars (|Cowie et al Jl99fj : lBarger ^BE^T nzalez et all 



120081 ). A new analysis in progress by Wall et al. (2009 in preparati on) uses a Bayesian 
modelling process similar to that described by IWall et al. and this pseudo- 



downsizing effect is very clear for FSRQ, as shown in Fig. 1101 

There remains further need for complete redshift information on faint samples of 
steep-spectrum radio sources to clarify their high-redshift evolution. An important 
step in this direction has been the Combined EIS-NVSS Survey Of Radio Sources 
(CENSORS. iBrookes et al1l200ol ) that included spectroscopic observations of 143 out 
of a total of 150 sources with S1.4GHZ > 3.8 mjy. Of these, 137 form a complete sample 
to a flux-density limit of 7.2 mjy. The result ing redshift distribution agrees well with 
the distribution in Fig. 28 of lCondonl (jl984al ) but is not well reproduced by any of the 
iDunlop fc Peacockl ( 1990l ) models (see, e.g., Fig. These data promise substantial 



improvement in this field. 
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N t01 = 18.79 deg" 




Fig. 11 Reds hift distribution of s ources brighter than 10 mjy at 1.4 GHz in the CEN- 
SORS sample (Brookes ct al. 2008), compared with that predicted by the PLE model by 
iDunlop &: Peacock! ||1990| . dashed line) and with the fit of eq. (I26C As the model applies to 
the AGN population, the two sources classified as starburst galaxies (nos. 95 and 124) were 
excluded from the histogram. The Kaplan-Meier estimator was used to take into account the 
lower limits into account. 



5.4 Unified evolutionary schemes 



lUrrv fc Padovanil l| 19951 ) carried out a thorough examination of the most widely ac- 
cepted version of the unification scheme encompassing steep- and flat-spectrum ra- 
dio sources. This scheme is based on the premise that relativistic beaming of lobe- 
dominated, steep-spectrum, moderate radio power FRI and high-radio-power FRII 
radio galaxies gives rise to core dominated, flat-spectrum BL Lac objects and radio 
loud quasars, respectively, when the line of sight is close to the jet axis (Fig. [TJ. Urry 
& Padovani based their analysis on a comparison of luminosity functions for the dif- 
ferent objects, and they showed that with reasonable assumptions for the beaming pa- 
rameters, these were in agreement. (We remind the reader that in the case of beamed 
emission, the true luminosities are lower than those inferred from the observed fluxes 
assuming isotropic emission by a factor oj/47r, where ui is the solid angle of the beaming 
cone(s). For example, in the case of twin beams each with semi-aperture angle of 7°, 
lu/Aty ~ 0.094.) 

More recent analyses ( Liu fc Zhandl2007l ; ICara fe Listed [2008). using observations 
of the jet k inematics and the apparent superluminal speeds for a complete sample 
(|Listerl(200sl ), have confirmed the general validity of the scheme and have improved the 
accuracy in parameter determination. The assumption that the observed luminosity 
function of FR11 galaxies has the same power-law shape as the intrinsic luminosity 
function of radio-loud quasars w as successfully tested using a maximum likelihood 
method. IPadovani fc Urrvl (jl992T l had originally developed their method of comparing 
the calculated beamed luminosity functions of flat-spectrum quasars and BL Lac ob- 
jects with the observed ones before their comprehe nsive review paper. N ow the observed 
distribution of Lorentz factors of relativistic jets ( Hovatta et al.ll2008l ) is found to be 
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in good agreement with their estimate. The observed distribution of viewing angles 
also tu rned out to be consistent with predictions of their unified model. iHovatta et al.l 
(2008|) found that the transition from blazars to or dinary radio loud quasar s occurs at 
a viewing angle of 15° to 20°, to be compared with lUrrv fc Padovanil (jl995h 's estimate 



of 14° 



Wall fc Jackson] (|l997h and I Jackson fc Walll l|l999h derived an evolutionary model 



aimed at explaining the behaviou r of both the flat- and s teep-spectrum populations 
within a unified scheme. Following lUrrv fc Padovanil (|l995h . they assumed BLLacs to 



be the beamed versions of FRI galaxies, and FSRQs to be the beamed versions of FRII 
galaxies. However the analysis differed from that of Urry & Padovani in deriving model 
parameters directly from the data rather than by comparing luminosity functions. As a 
first step, Jackson & Wall derived evolution models separately for FRI and FRII radio 
galaxies, using data (counts and redshifts) solely from low-frequency surveys in order 
to avoid beamed objects and to establish space-density behaviour for these parent 
populations. The evolution models they derived happened to have redshift cutoffs, 
and for the FRI galaxies the approximation of space density constant with epoch was 
adopted. These both gave good fits to the low-frequency data, but neither feature is 
essential to the outcome of the experiment. Together with observed ratios of beamed 
(core) flux to unbeamed flux, Jackson & Wall then used a Monte Carlo process to orient 
statistically-large samples of FRI and FRII sources randomly to the line of sight. Using 
a grid of beaming parameters (range of Lorentz factors, torus opening angles), they 
then calculated the number of beamed objects (BLLacs and FSRQs) produced at 
different flux densities. Identifying these objects with the cm-excess sources found in 
cm-wavelength surveys, they then closed the loop by using the higher-frequency source 
counts (primarily at 5 GHz) to find the permissible range of beaming parameters. 
The process was able to reproduce the exact form of the higher- frequency counts with 
their broader bump (§ 15. 2p as well as t orus opening angles and Lorentz factors in 



reasonable agreement with observations. iLiu fc Zhangl (2007) found that the Lorentz 
factor distribution is much steeper for low-redshift (z < 0.1) low-luminosity sources 
than for the more powerful, high-z ones, although the uncertainties are large. Therefore 
the most extreme relativistic jets are rarer in the low-z population. This indicates that 
the low- and high-redshift groups are likely t o be from different par ent populations, 
consistent with the dual-population scheme of Ijackson fc Wall (| 19991 ). 



The lJackson fc Walll jl999h unification scheme was adopted bv lWilman et al] (|2008l ) 



to reproduce the observed variety of radio-loud AGNs, including radio galaxies, steep- 
and flat-spectrum quasars, and strongly beamed sources such as BLLacs and FSRQs. 
These authors produced a simulation of a sky area of 20 x 20 deg 2 out to z = 20, 
and down to a flux density of lOnJy at several frequencies from 151 MHz to 18 GHz. 
The model comprises 5 source populations with different evolutionary properties: FRI 
and FRII AGNs, "radio-quiet" AGNs (defined as all X-ray - selected AGNs), and 
star-forming galaxies in quiescent and star-bursting phases. The simulation includes 
redshift-dependent size distributions both for radio-loud AGNs and for star-forming 
galaxies. Clustering is modeled by attributing to each population an effective halo 
mass and computing the corresponding bias factors, b(z). The simulations have been 
built with the SKA in mind and will serve to inform design of the SKA, design of 
analysis software and design of observing programmes. 
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5.5 Special classes of sources 

5. 5. 1 GHz peaked spectrum ( GPS) sources 

GPS sources are powerful (logLi.4GHz ~ 32 ergs -1 Hz -1 ), compact (;$ lkpc) radio 
sources with convex spectra peaking at GHz frequencies (see lO'Dealll99Sl . for a compre- 
hensive review). They are identified with both galaxies and quasars. It is now widely 
agreed that GPS sources correspond to the early stages of the evolution of powerful 
radio sources, when the radio-emitting region grows and expands within the interstel- 
lar medium of the host galaxy, befo re plunging out into the intergalactic medium and 
becoming an extended ra dio source dFanti et al.lll995l ; lReadhead et al.| [l996; Bcge lmanl 
Il996l ; ISnellen et ailboool ) . Conclusive evidence that these sources are young came from 
VLBI measurements of propagation velociti es. Speeds of up to ~ 0.4c were mea- 
sured, implying dyn amical ages ~ 10 years (|Polatidis et all 1 19991 ; iTavlor et ai1l2000l ; 



iTschager et alj 200Cf ). The identification and investigation of these sources is therefore 



a key element in t he study of the early e volution of radio-loud AGNs. 

The model by |Pe Zotti et al.l (|2000h implies that extreme GPS quasars, peaking 
at v > 20 GHz, should comprise a substantial fraction of bright radio sources in the 
WMAP survey at v ~ 20 GHz, while GPS galaxies with similar f poa k should be about 
10 times less numerous. For a maximum rest-frame peak frequency v p .i = 200 GHz, 
the model predicts about 10 GPS quasars with S30GHZ > 2 Jy peaking at > 30 GHz 
over the 10.4 sr at |6| > 10°. Although the number of quasars with spectral peaks at 



> 30 GHz in the WMAP survey is consistent with this, when additional data ( TrushkirJ 
2003) are taken into account such sources look more like bl azars caught durin g a flare 
that is optically thick up to high frequencies. Furthermore, iTinti et al.l (|2005r ) showed 
that most (perhaps two thirds) of the quasars in the sample of High F requency Peaker 
(HFP , GPS sources peaking above a few GHz) candidates selected bv lDallacasa et all 



(2000) are likely to be blazars, while all the 10 candidates classified as galaxies are 



consistent with being truly young sources. T his conclusion was stren gt hened by the 
VLBA varia bility and polarization studies of [Torniainen et al.l ( 20051 ) . lOrienti et al.l 
< 200d . l2007h . and lOrienti fc Dallacasal j200$ ). 



An implication of these results is that the samples of confirmed GPS quasars are too 
small to allow meaningful stu dy of their evolutionary properties. The situation is some- 
what better for GPS galaxies. ITinti fc de Zo tti ( 200(f) found that the observed redshift 
and peak frequency distributions of these sources can be satisf actorily accounted for i n 
the framework of the self-similar expansion model proposed bv lBegelmanl (| 19961 . 1 19991 ). 
According to this model, the properties of the sources are determined by the interaction 
of a compact, jet-driven, over-pressured, non-thermal radio lobe with a dense interstel- 
lar medium. Fits of the redshift and peak frequency distributions require a decrease of 
the emitted power and of the peak luminosity with source age or with decreasing peak 
fre quency, consistent w ith expectations from Begelman's model, but at variance with 
the lSnellen et al.1 {2000) model. 



5.5.2 Late stages of AGN evolution 

Late stages of the AGN evolution, characterized by low radiation/accretion efficiency, 
were brought into sharper focus by the discovery of ubiquitous, moderate-luminosity 
hard X-ray emission from nearby ellipticals. VLA studies at high radio frequencies (up 
to 43 GHz) have shown, albeit for a limited sample of objects, that all the observed 
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compact cores of ellip tical and SO galaxies have spectra rising up to ~ 20-30 GHz 
( di Matteo et al.lll999h . 

There is gro wing evidence that essent ially all massive ellipticals host super-massive 
black holes (e.g. iFerrarese fc Fordll2005l ). Yet nuclear activity is not observed at the 
level expected from Bondil (j 1952h spherical accretion theory, in the presence of ex- 
tensive hot gaseous h alos, and for the usually- assumed radiativ e efficiency ~ f0% 
( di Matteo et alJll999h . However as proposed bv lRees et al.l (|l982T ). the final stages of 
accretion in elliptical galaxies may occur via Advection- Dominated Accretion Flows 
(ADAFs), characterized by a very low radiative efficiency ( Fabian fc Reeslll995l) . The 
ADAF scenario implies strongly self-absorbed thermal cyclo-synchrotron emission due 
to a near-equipartition magnetic field in the inner parts of the accretion flows, most 
easily detected at cm to mm wavelengths. However the ADAF scenario is not the only 
possible explanation of the data, and is not problem-free. Chandra X-ray observations 
of Sgr A at the Galactic Center are suggestive of a co nsiderably lower accretion rate 
compared to Bondi's prediction (IBa ganoff et al.ll2003T ). so that the very low ADAF 
radiative efficiency may not be required. 

A stronger argument against a pure ADAF in the Gal actic Center is that the 
emiss ion is stron gly polarized at mm/s ub-m m wavelengths dAitken et al.ll2000l ; lAgoll 



l2000l l - Moreover Idi Matteo et al.l (jl999T ) and |Pi Matteo et al.l (|200lh found that the 
high-frequency nuclear radio emission of a number of nearby early-type galaxies is 
substantially below the prediction of standard ADAF models. The observations are 
more consistent with the adi abatic inflow-outflow solutions (ADIOS) developed by 
iBlandford fc Begelmanl (|l999T ). whereby the energy liberated by the accretion drives 
an outflow at the polar region. This outflow carries a considerable fraction of the mass, 
energy and angular momentum available in the accretion flow, thus suppressing the 
radio emission from the inner regions. Both the intensity and the peak of the radio 
emission depend on the mass loss rate. 

Tentative estimates of the counts due to these source s and o f the associated small- 
scale fluctuations were presented by IPierpaoli fc Pernal (|2004T l. As shown in Fig. [5] 
their model A assumes a local luminosity function (upper boundary of the cross- 
hatched area) higher than current estimates of the local luminosity function of all 
flat/inverted spectrum sources. Consistency is obtained only for their minimal model 
B (lower boundary of the cross-hatched area). 



6 Evolutionary models: star-forming galaxies 

6.1 Star- forming, normal and sub-mm galaxies 

The radio emission of star-forming galaxies correlates with their star forma tion rate, as 
demonstrated by the well - established tight correlation with far-IR emissi on dHelou et al.l 
ll985l ; lGavazzi et al.lll98d ; ICondorJl99alGarrettll2002r i. lYun etafl (|200lh found that the 
overall trend in the range Liy (60/im) ~ 10 30 -10 32 ' 5 ergs _1 Hz -1 is indistinguishable 
from a linear relation: 

ii.4GHz = 1.16 x 10~ 2 L„(60^m). (14) 

Galaxies with L v (&Qfim) < 10 30 ergs -1 Hz -1 are found to have radio to far-IR lu- 
minosity ratios systematically lower than those given by eq. (|14[) . The apparent de- 
viation from linearity in the radio/far-IR correlation at low luminosities is supported 
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by a comparison of 60 jj,m and 1.4 GHz local luminosity functi ons dYun et al 



Best et al.ll2005h. Simply shifting the 60 /im luminosity function ijSaunders et al 



2001 



1990 



Takeuchi et al.l [2003j) along the luminosit y axis according to eq. (1141) yields a g ood 
match to the radio luminosity function 1 Best et al.l 120051 ; iMauch fc Sadler! 120071 ) for 
^l.4GHz ~ 10 28 ergs~ 1 Hz^ 1 . At yet lower luminosities, however, the extrapolated lu- 
minosity function lies increasingly above the observed one. Full agreement is recovered 
(Fig. [U by replacing eq. fl4]) with 



^l.4GHz 



= 1.16 x 10~ 2 L h 



L;,(60/im) 



-3.1 



L u (60/j,m) 



(15) 



in which L b = 8.8 x 10 ZM ergs _i Hz~\ 

While a radio/far-IR correlation is expected since young stars are responsible both 
for dust heating and for the generation, via supernova explosions, of synchrotron emit- 
ting relativistic electrons, a clear explanation of its tightness and of its linearity over 
a large luminosity range is still missing. A decrease of the Z/! 4Qjj z /L„(60/im) ratio 
with increasing far-IR luminosity is expected from the increase of the effective dust 
temperature, T^, with luminosity (jBlain fc Longainll996h . For a galaxy like the Milky 
Way, the far-IR SED peaks at 170 n m, whereas for an Ul tra Luminous Infrared Galaxy 
(ULIRG) it peaks at about 60 /J.m (jLaeache et al.l 120051 ). This factor of 3 increase in 
temperature for a factor ~ 10 3 increase in luminosity corresponds to oc ipjj^- If the 
radio luminosity is proportional to the global far-IR luminosity, this increase in dust 
temperature results in a decrease of the Li .4GHz/^(60/xm) ratio by a factor of 2.5-3. 

On the other hand, there are different contributions to the global far-IR luminosity. 
In Luminous and Ultra Luminous Infrared galaxies, the emission is dominated by 
warmer dust, associated with star-formation, while infrared "cirrus" emission, heated 
by older stars, becomes increasingly important in galaxies with lower and lower star- 
formation rates. The latter component may be weakly correlated with radio emission, 
if at all. Moreover, in very low luminosity galaxies interstellar magnetic fields may be 
so weak as to let synchrotron emitting electrons escape into intergalactic space or to 
lose energy primarily via inverse Compton scattering of CMB photons. These processes 
may over-compensate the effect of decreasing dust temperatur^B 

Anyway, the tight empirical relationship between radio and far-IR luminosities 
for star-forming galaxies allows us to take advantage of the wealth of data at far- 
IR/sub-mm wavelengths to derive the radio evolution properties. We expect a different 
evolution for starburst and normal late-type galaxies as the starburst activity is likely 
triggered by interactions and mergers that were more frequent in the past, while in 
normal galaxies the star-formation rate has probably not changed m uch over their life- 
times. The bulk of the sub-mm counts measured by SCUBA surveys |Scott et al.ll2006l ; 
ICoppin et al]|2006h is due to yet another population, the sub-mm ga laxies (SMGs), 
proto -spheroidal galaxies in the process of forming most of their stars i Granato et al.l 
120041 ). 

There have been a number of attempts to model the evolution of star- forming galax- 
ies and in particular to account for the a pparent intrusion of this popula tion into the 
source counts at 5i.4GHz < 1 m Jyi see e.g. lKing fc Rowan- Robinson! 1 2004j). A straight- 



forwar d extrapolation to radio frequencies of the evolutionary models by Negrello et al.l 
(2007) for the three populations (normal, starburst and sub-mm galaxies), exploiting 



5 We are grateful to J. Condon for enlightening comments on this issue. 
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eq. (|15p and the SEDs of NGC 6946 for normal late-type galaxies and of Arp220 for 
starburst and proto-spheroidal galaxies, yields the curves shown in Figs. HJ{6] nicely 
reproducing the counts at tens of /iJy levels. We note however that new observational 
data, some of which is described in § 12.31 may permit substantial refinement of these 
models. 

The cross-over between synchrotron plus free-free emission prevailing at cm wave- 
lengths, and thermal dust emission, generally occurs at A ~ 2-3 mm (in the rest frame), 
so t hat at frequencies of tens of GHz there are contributions from both components 
(see |Pe Zotti et ail 120051) . 



6.2 Radio afterglows of 7-ray bursts (GRBs) 

The afterglow emission of GRBs can be modelled as synchrotron emission from a decel- 
erating blast wave in an ambient medium, plaus i bly the interste llar medium of the host 
galaxy (jWaxmanll 19971 ; IWiiers fc Galamalll999l ; lMeszaroslll999l '). The radio flux above 
the self-absorption break at ^ 5 GHz, is proportional to j/ 1 ' 3 up to a peak frequency 
that decreases with time. This implies that surveys at different frequencies probe dif- 
ferent phases of the expansion of the blast wave. Owing to their high brightnesses, 
GRB afterglows may be detected out to exceedingly high redshifts and are therefore 
important tracers of (a) the early star formation in the Universe, and of (b) the absorp- 
tion properties of the intergalactic medium across t he reionization p hase. Estimates of 
the counts of GRB afterglows have been made by ICiardi fc Loebl ( 2000t ) . who found 



that at a fixed time-lag after the GRB in the observer's frame, there is only a mild 
change in the observed flux density at radio wavelengths with increasing redshift. This 
stems in part from the fact that afterglows are brighter at earlier times and that a 
given observed time refers to an intrinsic ti me in the sour c e fram e that is earlier as 
the source redshift increases. According to ICiardi fc Loebl i|2000h estimates, a large 



area survey at ~ 1 cm to a flux limit ~ 1 mj y should discover some GRBs (see also 
ISeaton fc Partridgdl200ll ; be Zotti et "al"ll2005l ) . Predictions of Ciardi & Loeb's models 
at 1.4 and 3 GHz are shown in Figs. [S]and|51 



7 The radio background 

The radio background provides a key constraint on the counts of sources too faint to be 
individually detected. At frequencies below 1 GHz, the extragalactic radio background 
is swamped by the much more intense Galactic emission, primarily the synchrotron 
emission from diffused and inte grated super nova remnants. Estimates at meter wave- 



lengt hs made over 40 years ago (|Bridlelll967h using the T—T plot method (jTurtle et al 



Il962ft yielded an antenna temperature of the extragalactic background at 178 MHz of 
T"bkg = 30±7 K, about one third of the minimum total sky brightness at that frequency. 
Subtraction of the Galactic emission assuming that it scales as csc|6|, (b = Galactic 
latitude) is very inaccurate, as (i) Galactic emission towards the Galactic poles is not 
removed, and (ii) there are major features such as the North Galactic Spur which follow 
no such law. In fact the morpho logy is complex; see for example the superb map at 
408 MHz bv lHaslam etafl jl982j ). 



A n independent estimate of the background intensity was obtained bv lClark et al.l 
(1970), exploiting the low-frequency measurements obtained with the Radio Astron- 
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Fig. 12 Estimates of the extragalactic radio background at different frequencies. The black 
symbols refer to estimates exploiting the methods mentioned in the first two paragraphs of §[7] 
while the red symbols refer to estimates by Fixsen et al. (2009) either using ARCADE 2 data 
(squares) or re-analyzing published data from large area surveys (x). The dashed line shows 
the extragalactic background spectrum yielded by the models fitting the counts in Figs. I4"tf7l 



omy Explorer (RAE-1) satellite. While the T-T plot method exploits the isotropy of 
the background and its different spectral index to separate it from Galactic emission, 
low-frequency measurements exploit the strong attenuation of the extragalactic ra- 
diation below 1 MHz due to free-absorption by electrons in the interstellar medium. 
Measurements at these frequencies can therefore be used to single out the Galactic-disk 
component. Extrapolating to higher frequencies and subtracting from measurements 
of the total flux Clark et al. were able to obtain an estimate of the extragalactic back- 
ground intensity. Obviously the method works best in regions of low Galactic emission, 
such as the 'north halo minimum' region (/ ~ 150°, b ~ 50°). Once again the isotropic 
component was identified to be about one-third the minimum total brightness observed 
at 100 MHz. The spectral index appeared to be similar to the average spectral index 
observed for extra-galactic sources, suggesting that the isotropic component does rep- 
res ent the extragalactic background rather than an isotropic halo of the Galaxy. 

iGervasi et al.l (|200ct ) calculated the brightness of the isotropic background antici- 
pated from unresolved extragalactic source, using modern compilations of source counts 
and fitting smooth functions to these counts. Their results range from Tf, = 38600 K 
at 151 MHz to 0.41 K at 8.44 GHz; over this range the - frequency law is close to 
a power law. 

Interest in the radio background was recently revived by the results of the second- 
generation balloon-borne experiment ARCADE-2 (Absolute Radiometer for Cosmol- 
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ogy, Astrophysics, and Diffuse Emissi o n). A fter subtracting a model for the Galactic 
emission and the CMB, iFixsen et aD (^OO&h found excess radiation at 3 GHz about 



5 times brighter than the estimated contribution from extragalactic radio sources, as 
calculated by Gervasi et al. and in the present work (see Fig, I12[) . From a re-analysis of 
several large-are a surveys at l o wer fr equencies to separate the Galactic and extragalac- 
tic components, IFixsen et all l|2009h obtained an extragalactic background power-law 



spectrum of T = 1.26 ± 0.09 K (v/u Q y 2m±0M , with v Q = 1 GHz from 22 MHz to 10 
GHz, in addition to a CMB temperature of 2.725 ±0.001 K. These results are compared 
to earlier estimates in Fig. 1121 where we show the background brightness, I v , in Jy/sr, 
as a function of frequency. I v is related to the antenna temperature, T a , by 



c2l » onr ,„-6/ 



2ks^ VGHz/ Jy/sr 



-23 , 



As shown by the figure, the ARCADE results are inconsistent with earlier measure- 
m ents. We note that t he antenna temperature at 81.5 MHz implied by the power-law fit 
of IFixsen et al. I l|2009h . 854 K , exce eds the minimum total sky brightness temperature 



of 680 K measured by bridle! (|l967l ) 



8 Sunyaev-Zeldovich effects 



The Sunyaev & Zeldovich (SZ) effect jSunvaev fc Zeldovich|[l972T ) arises from the in- 



verse Compton scatter of CMB p h otons against hot electrons. For a comprehensive 
background review see iBirkinshawl (|l999l ). The CMB intensity change is given by 

V = (fec)2 W( - X) ( 7) 



where T CMB = 2.725 ± 0.002 K l|Mather et al.l 1 19991 ) is the CMB temperature and 
x = his/kT CMB . 

The spectral form of this "thermal effect" is described by the function 

g(x) = x 4 e x [x ■ coth(a;/2) - 4]/(e x - l) 2 , (18) 

which is negative (positive) at values of x smaller (larger) than xq = 3.83, corresponding 
to a critical frequency vq — 217 GHz. 
The Comptonization parameter is 

I kT e Ue fa Mq\ 
J m e c z 

where m e , n e and T e are the electron mass, density and temperature respectively, or is 
the Thomson cross section, and the integral is over a line of sight through the plasma. 

With respect to the incident radiation field, the change of the CMB intensity across 
a galaxy or a cluster can be viewed as a net flux emanating from the plasma cloud, 
given by the integral of intensity change over the solid angle subtended by the cloud 

AF V = I AI v dQ oc Y = / ydQ. (20) 
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In the case of hot gas trapped in the potential well due to an object of total mass M, 
the parameter Y in eq. (|20p . called integrated Y-flux, is proportional to the gas-mass- 
weighted electron temperature (T e ) and to the gas mass M g — f g M: 

Y(xf g (T e )M. (21) 

At frequencies below 217 GHz, the Y-flux is negative and can therefore be distinguished 
from the positive signals due to the other source populations. 



8.1 Sunyaev-Zeldovich (SZ) effects in galaxy clusters 

The SZ effect from the hot gas responsible for the X-ray emission of rich clusters of 
galaxies has been detected with high signal-to-noise and even imaged in many tens of 
objects dCarlstrom et al.ll2002l; iBenson et al.ll2004|; [ Jones et al.|[2005l ; iBonamente et ail 
2006; Hal verson et al.ll2009l ; IStaniszewski et al.l2009T ). Detailed predictions of the counts 
of SZ effects require several ingredients, generally not well known: the cluster mass func- 
tion, the gas fraction, the gas temperature and density profiles. All these quantities 
are evolving with cosmic time in a poorly-understood manner. Therefore, predictions 
are endowed with substantial uncertainties. Current models generally assume a self- 
similar evolution of the relationships betwee n the main cluster param eters (mass, gas 
temperature, gas fraction, X-ray lumin osity; Bonamente et alJ 20 0^). Several predic- 
tions for the SZ counts are available (e. g. de Luca et al .11 19951 ; I IColafrancesco et al ] |l997l ; 
|Pe Zotti et~aHl2005l ; lchamballu et alJbOQSh . SZ~ maps have be en constructed, mostly 
using the o u tput o f hydrodynamical cosmological sim ulations, bv lGeisbiisch et all (2005 ) , 
IPace et"all l|2008l ). IWaizmann fc Bartelmanr] (2J309), amongst others. 

Our understanding of the physics of the intra-cluster plasma is expected to improve 
drastically than ks to ongoing and forth coming SZ surveys such as those with the South 
Pole Telescope dCarlstro m et al ]|2009h . the Atacama Cosmolog y Teles cope dKosowskv 

2006 ) . APEX-S zlDobbs et alJbood). AMI (|Zwart et alj|2008r i . SZA jMuchovei et al. 

2007 ) . OCRA_p |Lancaster et alfeOQTT ). and the Planck mission (|The Planck Collaboration! 



2006) 



8.2 Galaxy-scale Sunyaev-Zeldovich effects 

The formation and early evolution of massive galaxies is thought to involve the release 
of large amounts of energy that may be stored in a high- pressure proto-galactic plasma, 
producing a detectable SZ effect. (Note that the amplitude of the ef fect is a measure 
of the pressure of the pla sma.) According to the standard scenario i Rees fc Ostrikerl 



1977l ; IWhite fc Reesll97Sr ). the collapsing proto-galactic gas is shock- heated to the virial 
temperature, at least in the case of large halo masses ( > 10 12 M^. lDekel fc Birnboiml 
2006). Further important contributions to the gas thermal energy may be produced 
by supernova explosions, winds from massive young stars, and mechanical energy re- 
leased by central super-massive black-holes. The corresponding SZ signals are poten- 
tially a direct probe of the processes governing the early phases of galaxy evolution 
and on the history of the baryon content of gala xies . They have been investigated 
unde r a variety of assumptions by many authors (fohlll999l: Nataraian & Sigurdsson 
19991 ; lYamada et alJll99alAghanim et alJbood : iMaiumdar et al.ll200ll ; iPlatania et al.l 
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20021: lOh et al]|2003l: lLapi et al.ll20ol iRosa-Gonzalez et al.ll2004l: foe Zotti et al.ll2004l : 
Chatterie^fcKosowskvl 2007 : iMassarchet'ld* J l2008bl : IChatteriee etaLl 2008). 



Widely different formation modes for present day giant spheroidal galaxies are being 
discussed in the literature, in the general framework of the standard hierarchical clus- 
tering scenario. One mode ( Granato et aljliool ; lLapi et al.|[200r3 ; ICook et al.ll2o"o9h has 



it that these galaxies generated most of their stars during an early, fast collapse featur- 
ing a few violent, gas rich, major mergers; only a minor mass fraction may have been 
added later by minor mergers. Alternatively, spheroidal galaxies may have acquired 
most of their stars throu gh a sequence of, mostly dry, mergers 1 De Lucia fc Blaizotl 
120071 : iGuo fc Whitell200Sl ). 

The second scenario obviously predicts far less co nspicuous galaxy-scale S Z signals 
that the first one. In the framework of the first scenario iMassardi et~aLl (|2008bh find that 



the detection of substantial numbers of galaxy-scale thermal SZ signals is achievable 
by blind surveys with next generation radio telescope arrays such as EVLA, ALMA 
and SKA. This population is detectable even with a 10% SKA, and wide-field-of-view 
options at high frequencies on any of these arrays would greatly increase survey speed. 
An analysis of confusion effects and contamination by radio and dust emissions shows 
that the optimal frequency range is 10-35 GHz. Note that the baryon to dark matter 
mass ratio at virialization is expected to have the cosmic value, i.e. to be about an 
order of magnitude higher than in present day galaxies. Measurements of the SZ effect 
will provide a direct test of this as yet unproven assumption, and will constrain the 
epoch when most of the initial baryons are swept out of the galaxies. 



9 Wide area surveys and large scale structure 

Extragalactic radio sources are well suited to probe the large-scale structure of the 
Universe: detectable over large cosmological distances, they are unaffected by dust 
extinction, and can thus provide an unbiased sampling of volumes larger than those 
usually probed by optical surveys. On the other hand, the ir 3D space-distribution can 
be recovered only in the very local Universe (z ;$ 0.1; see lPeacock fc NicholsorJll99ll ; 

iMaeliocchetti et al.ll20o3 l because the majority of radio galaxies detected in the avail- 
able large-area, yet relatively deep, surveys, carried out at frequencies < 1.4 GHz, have 
very faint optical counterparts, so that redshift measurements are difficult. As a result, 
only the angular (2D) clustering can be measured for the entire radio AGN popula - 
tion. High-frequency surveys have much higher identification rates |Sadler et alJbooih . 
suggesting that this difficulty may be overcome when such surveys cover sufficient sky 
and are linked to wide-area redshift surveys. 



9.1 The angular correlation function and its implications 

Just the basic detection of clustering in the 2D distribution of radio sources proved to be 
extremely difficult (e.g. IWebsterlll976l ; ISeldner fc Peeblej[l98ll : IShaver fc Pierrell 19891 ) 
because at any flux-density limit, the broad luminosity function translates into a broad 
redshift distribution, strongly diluting the spatial correlations when projected onto the 
sky. Only with the advent of d eep radio surveys cove ring larg e areas of the sky, F IRST 
dBecker et al.ll 19951). WENS S jRengelink et all 19971 ). NVSS (|Condon et al.lll99Sl ). and 
SUMSS jMauch et all [20031 ) . did it become possible to detect the angular clustering 
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of t hese objects with high statistica l significance: s ee Cress fc Kamion cowski|j[l998h 
and IMagliocchetti et al.1 Jl998l . Il999 ) : iBlake fc Walll d2002al) for FI RST; iBlake fc Wal 
J2002bllai) a n d lOverzier et al.l (|2003h for NVSS; iRengelinkl l| 19991 ) for WENSS; and 
IBlake et all (|2004bh . Even then there remained difficulties of interpretation due to 
spurious correlation at sm all angular scales cau sed by the multiple-component nature 
of extended radio sources <|Blake fc Wallll2002bh : the raw catalogues constructed from 
these large surveys list components of sources rather than single 'assembled' sources. 
Amongst the cited surveys, NVSS is characterized by the most extensive sky coverage 
and can thus provide the best clustering statistics, despite its somewhat higher com- 
pleteness limit (~ 3mJy vs ~ 1 mjy of FIRST). The two-point angular correlation 
function w(6), measured for NVSS sources brighter than 10 mjy, is well described by a 
power-law of slope —0 .8 extending from ~ 0.1 degrees up to scales of almost 10 degrees 
(|Blake fc Wai] l2002ah . A signal of comparable amplitude and shape was detected in 
the FIRST survey at the same flux density limit, on scales of up to 2-3 degrees (see 
e.g. IMagliocchetti et al.|[l99sl . Il999l) . while at larger angular separations any positive 
clustering signal - if present - is hidden by the Poisson noise. 

Most of the analyses performed so far wit h the aim of reproduci ng the clustering 

et alll2003h assumed a two- 
(7*/ Vn)~ 7 . The power-law shape 



es performed so tar wit n tne 
i.g. blake fc Waiill2002bfl ; [Ch 



of radio galaxies (see e.g. IrJlake fe Walill^UU2bl lal; lUverz u r 
point spatial correlation function of the form £rgM 

is in fact preserved when projected onto the sky l|Limberlll953l ). so that the observed 
behaviour of the angular correlation is well recovered. The correlation length ro was 
found to lie in the range 5-15 Mpc, the large range reflecting the uncertainties in both 
the redshift distribution of the sources and the time-evolution of clustering. Despite 
the wide range in measurement of ro, the above results suggest that radio galaxies are 
more strongly clustered than optically-selected galaxies. 

A deeper examination of the power-law behaviour of the angular two-point corre- 
lation function up to scales of the order of ~ 10° highlights interesting issues. Within 
the Cold Dark Matter paradigm of structure formation, the spatial correlation function 
of matter displays a sharp cut-off near a comoving radius of r ~ 100 Mpc, which at 
the average redshift for radio sources < z >~ 1, corresponds to angular separations 
of only a few (~ 1° — 2°) degrees. This is in clear contrast to the observations of 
the angular two-point correlation function. The question is how to reconcile the clus- 
tering properties of these sources with the standard scenario of structure formation. 
Some authors have tried to explain the large-scale positive tail of the angular correla- 
tion function w(9) as due to a high-density local popu lation of star-forming galaxies 
(|Blake et al.ll2004ah . Others (|Magliocchetti et al.lll999l ) suggested that the results can 
be reproduced by a suitable choice of the time-evolution of the bias parameter, i.e. 
the way radio galaxies trace the underlying mass distribution. The first hypothesis can 
be discarded on the basis of more recent determinations of the space density of local 
star- forming galaxies with a 1.4-GHz rad io counterpart (e.g. IMagliocchetti et al,| [2002; 
ISadler et al.il2002l ; iMauch fc Sadlerll2007h . Even the second approach, although promis- 
ing, suffers a number of limitations due to both theoretical modelling and quality of 
data then available. 

Theoretical predi ctions for the angular two-point correlation function of a given 
class of objects using iLimberl (|l953h equation 
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require two basic ingredients: their redshift distribution, N(z), i.e. the number of ob- 
jects brighter than the flux limit of the survey as a function of redshift, and the value 
of the bias factor as a function of redshift, b(z). In equation (122f) . r(Sz,8) represents 
the comoving spatial distance between two objects located at redshifts z and z + Sz 
and separated by an angle 8 on the sky. For a flat universe and in the small-angle 
approximation (still reasonably accurate for scales of interest here, 0.3° ^ 8 ;$ 10°) 



2 ( c 



Sz 



(dc(z)ey 



with 



e(z) = [n M (i + zf + n A ] 1/2 , d c ( z ) 



dz' 



(23) 



(24) 



On sufficiently large scales, where the clustering signal is produced by galaxies residing 
in distinct dark-matter halos and under the assumption of a one-to-one correspondence 
between sources and their host halos, the spatial two-point correlation function can be 
written as the product of the correlation function of dark matter, £pM, t imes the square 
of the bias parameter, b i Matarrese et etl.ll 19971 : iMoscardini et al. ll99Sh : 



£(r,z) = b (M eS ,z)£r> M (r,z). 



(25) 



Here, M c g is the effective mass of the dark matter haloes in which the sources reside 
and b is derived in the extended Press fc Schechteil ( 19741 ) formalism according to the 

pre scriptions of [Sheth fc Tormen ( 19991 ). 

iNeerello et alJ l|2006l ) adopted the Af(z) from bunlop fc Peacock] <|l99dl Vs pure lu- 
minosity evolution model. If the effective mass of the dark matter haloes in which the 
sources reside does not depend on cosmic time, as found for optically-selected quasars 
( Porciani etafl 120041 : [c loom et al.ll2004l ). ;he predicted angular correlation function 
badly fails to reproduce the observed one. This is because contributions to w(8) on 
a given angular scale come from both local, relatively close pairs of sources and from 
high-redshift, more distant ones. But for z ^ 1 angular scales 8^2° correspond to lin- 
ear scales where the correlation function is negative. Since the contribution of distant 
objects is overwhelming, we expect negative values of w(8), while observations give us 
positive values. 

The only way out appears to be a damping down of the contribution t o w(8 ) 
of high-z sources, and this can only be achieved through b(z). iNeerello et alJ (2006) 
found that the w(8) data can be reproduced by assuming an epoch-dependent effective 
mass proportional to the mass scale at which the matter-density fluctuations collapse 
to form bound structures. Such mass decreases with increasing redshift, thus abating 
the negative high- z contributions to w(8). This assumption may be justified - locally, 
AGN-powered radio galaxies are found mainly in very dense environments such as 
groups or clusters of galaxies, and the characteristic mass of virialized systems indeed 
decreases with increasing redshift. The best fit to the data was obtained for a high 
value of the local effective mass, M eS (z = 0) ~ 10 15 M©/h. However the CENSORS 



data I Brookes et alJlioOSl ) have shown that the reds hift distribution peaks at lower 
redshifts than predicted by iDunlop fc Peacock! |l990l ) PLE model (Fig. [Til)- Using a 
smooth description of the CENSORS redshift distribution 



N{z) = 1.29 + 32.37z - 32.89z 2 + 11. 13z 3 - 1.25z 4 , 



(26) 
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Fig. 13 Two-p oint angular correlatio n function of NVSS sources wi th 5i.4GHz > l OmJy 
as measured by Blake & Wall (2002b) compared with the model by [Ncgrcllo ct al. (2006, 
red curves) and wit h the updated model (blue curves) fitting the redshift distribution by 
IBrookes et all |2008). The dashed curves include the contribution of a constant offset e = 
0.0001 to w(9) in order to account for the effect of possible spurious density gradients in the 
survey. 



the best fit is obtained with a somewhat lower value for the local effective mass, 
M eS (z = 0) ~ 10 14 - 5 M Q /h (Fig. 



9.2 Integrated Sachs-Wolfe (ISW) effect 

The ISW effect describes the influence of the evolution of the gravitational potential 
in time- variable, linear, metric perturbations on CMB photons that traverse them. 
When the CMB photons enter an overdensity they are gravitationally blue-shifted, and 
they are red-shifted when they emerge. In an Einstein-de Sitter universe the density 
contrast grows as the linear scale, so that the gravitational potential associated with 
the mass fluctuation is independent of time. Hence the red- and the blue-shift exactly 
compensate each other and the net effect is zero. However, a non-zero effect arises if 
the gravitational potential decays, as in the case of an open universe when the effect of 
the space curvature is important, or when the dynamics of the universe are dominated 

by dark energy. 

As first pointed out bv lCrittenden fc Turokl ( 1996h . a promising way of probing the 



ISW effect is through correlating fluctuations in the Cosmic Microwave Background 
(CMB) with large-scale structure. Since the timescale for the decay of the potential 
is of the order of the present-day Hubble time, the effect is largely canceled on small 
scales, because photons travel through multiple density peaks and troughs. This is why 
surveys covering large areas of the sky and probing the large scale distribution up to 
z ~ 1 are necessary. 

A hi gh quality all-sky CMB temperature map has been provided by the WMAP 
satellite ( Bennett et al. I l2003l ; iHinshaw et "al1l2007l ; Iffinshaw fc Naevefeoosf ) . A particu- 



larly well-suited probe of the large-scale structure is the NVSS survey, and indeed this 



45 



has been extensively exploited to look for the ISW signal feouehn fc Crittenden! 2004 , 

20051 : [Pietrobon et "alJgOOdlMcEwen et alll2007l . l200Sl : IHo et all200Sl : lGiannantonio et al.1 



20081 ; iRaccanelli et al.ll200S 



The comparison of the correlations inferred from the data with model predictions 
requires once again the redshift distribution and the bias parameter as a function of 
redshift. All analyses carried out so far have used redshift distributions inconsistent 
with the CENSORS results. The product of the latter redshift distribution with the 
redshift-dependent bias factor best fitting the observed w(8) (see the previous sub- 
section), whose integral determines the amplitude of the ISW effect, peaks at redshifts 
where the contribution to the ISW signal in a /1CDM cosmology also peaks, namely 
z ~ 0.4. This means that the NVSS sample is very well suited to test the effects of 
dark energy on the growth of structure. The predicted cross-correlation power spectrum 
between the surface density fluctuations of NVSS sources and the CMB fluctuations 
expected for the "concordance" /1CDM cosmology turns out to be in good agreement 
with the empirical determination using the CMB map obtained from WMAP data. 
This conclusion is at odds with that of IHo et al. I {200I) who found that the WMAP 



3-year model predicts an ISW amplitude about 2<r below their estimate. Hence we 
suggest that the amplitude of the ISW cross-corre lation does not supp ort the case for 
new gravitational physics on cos mological scales llAfshordi et al"] |2008) or for a large 
local primordial non- Gaussianity (jAfshordi fc Tollevll2008I l. 



10 The future 

There are prospects for dramatic steps forward in radio and millimeter-wave astron- 
omy within the decade, thanks to a new generation of large to gigantic interferometers 
as well as refurbishment of old interferometers. Interferometric observations gain over 
single-dish observations not just through resolution but through improved sensitiv- 
ity, because correlation of the signals from the antennas can distinguish signal from 
noise and background. Long integrations become possible without the limitation of 
systematic errors. However, observing with interferometers requires careful set-up of 
the antenna array (and its parameters in software) to image sources correctly by mea- 
suring their flux densities on the appropriate angular scales. There is also the dreaded 
problem of 'missing flux' from (lack of) low-order harmonics in the spatial transform, 
corresponding to structure on the larger scales (§0)- In addition, the amount of post- 
processing required is large in comparison to single-dish measurements, to correct for 
the many instrumental and atmospheric issues, to convert the Fourier components into 
brightness images of radio sources, and (because of incomplete sampling in the Fourier 
plane) to apply algorithms to maximize image fidelity and dynamic range. Interferom- 
eters offer an additional advantage: on the increasingly noise-polluted surface of our 
planet, the processing can excise radio-frequency interference (RFI), which, even at the 
remote sites of future large arrays, would otherwise seriously compromise observations. 

The small beam size at millimeter wavelengths makes large-area deep surveys ex- 
tremely difficult because of the time penalty. New scanning techniques will need to 
be developed to perform such surveys. The next generation of interferometers, thanks 
to larger collecting areas, broader band widths and faster scanning capabilities are ex- 
pected to produce deeper surveys of large sky regions, both in total intensity and in 
polarization. Receiver advances have resulted in huge bandwidths (> 10 per cent) and 
very low equivalent noise temperatures over these bandwidths. To realize these gains in 
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sensitivity with interferometers requires development in correlator speed and in process- 
ing power. Moreover instead of single-pixel feeds, the development of focal plane arrays 
(FPAs, or phased-array feeds) looks to realize the long-standing dream of making near- 
full u se of the information brought to the focal plane (e.g. APERTIF. IVerheiien et ah! 
2008). There are further implications for correlator- and processor-power requirements. 

The focus of current effort in the radio-astronomy community is towards the re- 
alization of the Square Kilometer Array (SKA), the largest and most sensitive radio 
telescope ever. The SKA stands to be one of the iconic scientific instruments of the 
21st century. It will consist of an array of thousands of dishes, each 10-15 m in di- 
ameter, as well as a complementary aperture array - a large number of small, fixed 
antenna elements plus receiver chains arranged in a regular or random pattern on the 
ground. Between these two technologies, a frequency range of 100 MHz to 25 GHz will 
be covered. The collecting area will add up to approximately one million square me- 
ters, with baselines ranging from ~15 m to more than 3000 km. The SKA will require 
super-fast data transport networks (of order the current total internet capacity) and 
computing power far beyond current capabilities. Indeed the concept is only feasible 
if Moore's Law (the packing density of processing elements approximately doubling 
every two years) continues to hold; and this in itself requires revolutions in processor 
technologies. Site-testing has narrowed the choice to remote regions either in Western 
Australia or in South Africa. The telescope is expected to be fully operational after 
2020, but a 10% SKA may be operating as early as 2015. Many different technological 
solutions will be selected and integrated into the final instrument: they will represent 
the results of developing the so-called SKA pathfinders (ASKAP, MeerKAT, ATA, 
LOFAR, e-MERLIN, EVLA; www-page descriptions are readily available for each). 

These pathfinders carry the shorter-term excitement; they themselves represent 
leaps forward in observational capability. Most of them will be operational by 2015. 
They will realize dramatic improvements in survey speed and se nsitivity; ASKAP for 
example, is expected to produce a survey similar to the NVSS ( Condon et al.l ll998) 
in sky covera ge, but bettering it b y a factor of 50 in sensitivity and of 5 in angu- 
lar resolution ( Johnston et al.ll2008f ). These new tools will allow us to distinguish the 
star-forming population from the AGN population at low flux-density levels and to 
investigate source populations at extreme redshifts. 

The SKA itself will impact every area of astronomy and cosmology, from detection 
and mapping of planetary systems, study of individual stars, star clusters, pulsars, 
the structure of our Galaxy in both baryons and in magnetic field, through to normal 
galaxies, AGNs, pr oto-galaxies, and large-sc ale structure of the Universe. The com- 
pelling science (e.g. ICarilli fc Rawhngsll2004l, and updates on the SKA website) to be 
realized mak es irresistible reading. 

LOFAR ( Rottgering et af1l2006h . will open up the frequency window at the low end 
of the radio spectrum, below 240 MHz. LOFAR will survey the sky to unprecedented 
depths at low frequencies and will therefore be sensitive to the relatively rare radio 
sources that have very steep spectra, extreme luminosities and redshifts (§[T]). A unique 
area of investigation will be the search for redshifted 21cm line emission from the epoch 
of reionisation. 

At the other end of the radio band, the Atacama Large Millimeter Array (ALMA) 
is eagerly anticipated by the mm continuum and molecular-line community. Being 
built on a high (5000m), dry plain in the Atacama desert, northern Chile, it is an 
international project, a giant array of 50 12-m submillimetre quality antennas, with 
baselines of several kilometres. An additional compact array (ACA) of 12 7-m and 4 12- 
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m antennas is also foreseen. ALMA will be equipped with mm and sub-mm receivers 
covering ultimately all the atmospheric windows at 5000m altitude in ten spectral 
bands, from 31 to 950 GHz. The array will be operational by 2012 with a subset 
of the high-priority receivers. The steep rise of the dust emission spectrum at mm 
and sub-mm wavelengths implies that the K-correction compensates, at z ^ 0.1, for 
the dimming due to increasing distance (jBlain fc Longair|[l993T ). making the observed 



mm flux of dusty galaxies of given bolometric luminosity only weakly dependent on 
redshift up to z ~ 10. This makes ALMA the ideal instrument for investigating the 
origins of galaxies in the early universe, with confusion made negligible by the high 
spatial resolution. Using far-IR emission lines and CO rotational emission, ALMA will 
reveal the astrophysics of early phases of galaxy formation and provide the redshift of 
large numbers of obscured star-forming galaxies up to very large distances. This will 
enable us to establish the star-forming history of the universe, without the uncertainties 
caused by dust extinction in optical studies. 

Technological advances have resulted in upgrades of existing telescopes that revo- 
lutionize performance. For examples: 

1. The Australia Telescope Compact Array (ATCA), a 6 22m-dish array, has recently 
completed the upgrade of 7 mm receivers (working in the frequency range 30-50 
GHz), and the increase of the bandwidth from 2 x 128 MHz to 4 GHz (thanks to 
the new CABB system). These new capabilities together with its fast scan rate 
(15°/min at the meridian) will allow the extension of the Australia Telescope 20 
GHz (AT20G) Survey to higher frequencies or to lower flux densities. 

2. The Expanded Very Large Array (EVLA) is an upgrade of the sensitivity and 
frequency coverage of the VLA. When completed, it will use the 27 25m dishes of 
the VLA with 8-GHz bandwidth per polarization in the frequency bands 18-26.5, 
26.5-40, and 40-50 GHz. This is a 10 to 100-fold increase in sensitivity over the 
standard VLA. First observations will be in 2010; after full commissioning (2013) 
the (E)VLA is destined to remain at the forefront of radio astronomy for at least 
a decade. 

Finally, several new survey instruments for the Sunyaev-Zeldovich effect have either 
started operations or will shortly do so (see Ei l8.1j l. 



11 Conclusions 

For decades radio surveys have been a leading agent for extragalactic research, as 
testified by the breakthroughs they triggered, from the discovery of cosmic evolution, 
to quasars and the first high-z galaxies. They continue at the forefront of astrophysics 
and cosmology; e.g., via large-scale structure studies they pose challenges for buildup 
of the cosmic web; and via downsizing and AGN feedback deemed to produce this 
downsizing, they have come to the fore in modelling galaxy formation and evolution. 

But our physical understanding of the origin and evolution of the AGN-powered 
radio emission is still poor. While physical models for the cosmological evolution of 
galaxies and radio-quiet quasars have been progressing rapidly in recent years, the 
main progress on the radio side has been towards a phenomenological description of 
evolution of various radio AGN types. 

Even on the phenomenological side, there are aspects that are not fully understood. 
The epoch-dependent luminosity functions of galaxies and radio-quiet quasars are now 
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quite accurately determined up to high redshifts and there are attempts to provide 
physical explanations for the evidence for earlier formation of the more massive objects. 
Direct evidence for a substantial decline of the space densities of radio AGNs at z }t 2 
remains somewhat controversial, although modern evolutionary models accounting for 
the observed counts and redshift distributions do include such a decline. 

The origin of these uncertainties remains as it has been for the last 30 years - a 
lack of concerted effort to obtain complete redshift sets for radio AGN samples. The 
new redshift surveys (2dF, SDSS) have helped greatly in defining local space densities. 
But beyond two or three of the brightest samples such as 3CRR, there are no samples 
with complete redshift data. Even samples of 100 to 200 objects would suffice for most 
purposes, and could be easily obtained with 8- to lOm-class telescopes. Our deficiency 
in this regard is hig hlighted by the reliance of most analyses needing complete redshift 
information on the iDunlop fc Peacock! (|l990h model distributions - even 20 years on. 
Doubtless the planned new de ep optical and infrared wide-field surve ys such as PAN- 
STAR RS (|Hodapp et al.ll2004h and those with the VST and VISTA jArnaboldi et all 
2007) will help. However there will remain a need for individual pursuit of the faintest 
members in samples via deep imaging in different optical and IR bands and via fast 
spectrographs to complete the redshift information in samples of limited size. 

Sixty years since their discovery, radio AGNs remain at the forefront of astrophysics 
and cosmology. Our continued attempts to solve the mysteries which still surround 
them will doubtless lead to fresh discoveries of impact as great as those which have 
distinguished the first 60 years. 
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Appendix 

We present tabulations of the observed source counts, including when relevant, an estimate of 
the clustering contribution to the errors. The columns 'norm, cts.' and '+,-' give, respectively, 
S 2,5 dN/dS[Jy 1 ' 5 /sr], and its positive and negative errors. The column 'sampl. err.' gives, for 
surveys covering less than 25 deg 2 the contribution to such errors due to the sampling variance 
[eq. JB)]; such contributions are negligible for larger areas. 



Table 1: Euclidean normalized differential source counts at 150 MHz. 



log(S[JyJ) 


norm. cts. 


+ 




ref. 


-0.978 


765.720 


47.160, 


55.080 


mc90 


-0.908 


957.960 


59.400, 


68.760 


mc90 


-0.859 


1079.640 


66.600, 


77.760 


mc90 


-0.810 


1104.120 


68.040, 


79.560 


mc90 


-0.772 


1301.400 


80.280, 


75.600 


mc90 


-0.728 


1330.560 


82.080, 


86.400 


mc90 


-0.685 


1510.920 


93.240, 


98.280 


mc90 


-0.677 


1305.360 


90.720, 


88.920 


ha88 


-0.636 


1580.040 


97.561, 


91.800 


mc90 


-0.634 


1459.440 


89.640, 


68.040 


ha88 


-0.596 


1584.000 


84.960, 


50.760 


ha88 


-0.587 


1652.400 


128.161, 


107.639 


mc90 


-0.552 


1572.120 


84.240, 


50.401 


ha88 
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Table 1: continued 



gfSfJvl) 


norm 


cts. 


+ 


1 


ref. 


-0.533 


1947 


600 


135.360, 126.720 


mc90 


-0.515 


1879 


920 


85.680, 59.400 


ha88 


-0.471 


1922 


040 


118.080, 32.400 


ha88 


-0.467 


1918 


080 


133.560, 83.879 


mc90 


-0.428 


2070 


720 


94.680, 


79.560 


ha88 


-0.391 


2548 


080 


197.641, 


129.959 


mc90 


-0.385 


2230 


560 


101.880, 85.320 


ha88 


-0.347 


2569 


679 


117.360, 


115.560 


ha88 


-0.309 


2768 


040 


126.720, 


123.479 


ha88 


-0.304 


2472 


479 


172.081, 


160.919 


mc90 


-0.266 


2872 


800 


131.400, 


127.800 


ha88 


-0.228 


2937 


599 


157.321, 


109.799 


ha88 


-0.196 


3187 


080 


171.001, 


207.360 


mc90 


-0.184 


3117 


600 


191.519, 93.240 


ha88 


-0.141 


2981 


160 


206.639, 


153.720 


ha88 


-0.109 


3538 


079 


216.721, 


180.720 


ha88 


-0.071 


3618 


000 


251.999, 


184.679 


mc90 


-0.060 


4046 


400 


248.400, 


176.401 


ha88 


-0.016 


3726 


000 


230.401, 


189.000 


ha88 


0.027 


3866 


400 


208.800, 


194.399 


ha88 


0.076 


4136 


401 


284.399, 


237.602 


ha88 


0.130 


3726 


000 


259.199, 


242.640 


ha88 


0.190 


4388 


399 


338.401, 


251.998 


ha88 


0.190 


3837 


600 


237.600, 


223.201 


mc90 


0.255 


3895 


199 


331.201, 


280.800 


ha88 


0.321 


3808 


800 


356.400, 


274.681 


ha88 


0.391 


4978 


800 


464.401, 


388.799 


ha88 


0.467 


3614 


399 


424.801, 


336.238 


ha88 


0.565 


4068 


000 


450.001, 


457.201 


ha88 


0.668 


3805 


201 


482.400, 


454.321 


ha88 


0.783 


2823 


120 


502.919, 


500.760 


ha88 


0.919 


3349 


800 


595.799, 


613.080 


ha88 


1.147 


2927 


520 


626.041, 


554.399 


ha88 



Reference codes . ha 88: lHales et alj <| 19881 ^ : mc90: 
McGilchrist ct al. (1990|). 



Table 2: Euclidean normalized differential source counts at 325 MHz. 



log(S[Jyj) 


norm. cts. 


+ 




sampl. err. 


ref. 


-3.372 


25.400 


3.700, 


3.700 


0.044 


ow09 


-3.270 


25.900 


3.700, 


3.700 


0.044 


ow09 


-3.125 


16.900 


1.800, 


1.800 


0.029 


ow09 


-2.949 


17.200 


2.400, 


2.400 


0.030 


ow09 


-2.776 


19.300 


2.000, 


2.000 


0.033 


ow09 


-2.602 


27.300 


2.900, 


2.900 


0.047 


ow09 


-2.426 


30.800 


4.200, 


4.200 


0.053 


ow09 


-2.250 


58.000 


7.901, 


7.901 


0.100 


ow09 


-2.175 


33.935 


7.911, 


7.911 


0.241 


0088 


-2.077 


48.700 


9.400, 


9.400 


0.084 


ow09 


-1.873 


75.660 


9.784, 


9.784 


0.538 


0088 


-1.824 


89.000 


14.701, 


14.701 


0.153 


ow09 


-1.573 


166.577 


19.838, 


19.838 


1.185 


0088 


-1.523 


243.600 


43.702, 


43.702 


0.418 


ow09 


-1.272 


450.850 


51.078, 


51.078 


3.207 


0088 
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Table 2: continued 



log(5[Jy]) norm. cts. + , - sampl. err. rcf. 



-1.222 


495.600 


105.103, 105.104 


0.851 


ow09 


-0.971 


468.133 


86.201, 86.201 


3.330 


0088 


-0.921 


398.100 


162.501, 162.502 


0.684 


ow09 


-0.670 


1005.126 


213.841, 213.841 


7.150 


0088 


-0.218 


1535.170 


396.048, 396.048 


10.921 


0088 



Refere nce codes. ow09: 1 Owen et al.l Ull); 0088: Oort ct al. 
dl98SD. 



Table 3: Euclidean normalized differential source counts at 408 MHz. 



log(5[Jy]) 


norm. cts. 


+ 


_ 


sampl. err. 


ref. 


-1.857 


86.625 


22.647, 


22.647 


2.577 


be82 


-1.673 


124.875 


11.847, 


11.847 


3.715 


be82 


-1.488 


200.250 


23.275, 


23.275 


5.957 


be82 


-1.304 


304.875 


34.947, 


34.947 


9.069 


be82 


-1.120 


430.875 


46.790, 


46.790 


12.817 


be82 


-1.111 


437.400 


36.010, 


36.010 


0.835 


gr88 


-1.034 


441.000 


27.013, 


27.013 


0.842 


gr88 


-1.021 


491.908 


53.754, 


59.034 


_ 


ro77 


-0.936 


478.125 


58.020, 


58.020 


14.223 


be82 


-0.910 


619.861 


59.305, 


47.119 


_ 


ro77 


-0.903 


540.000 


19.827, 


19.827 


1.031 


gr88 


-0.757 


732.600 


23.442, 


23.442 


1.398 


gr88 


-0.752 


835.875 


115.215, 


115.215 


24.865 


be82 


-0.746 


762.302 


88.416, 


66.481 


_ 


ro77 


-0.604 


1039.737 


149.211, 


102.227 


_ 


ro77 


-0.602 


892.800 


27.054, 


27.054 


1.704 


gr88 


-0.568 


789.750 


148.125, 


148.125 


23.493 


be82 


-0.456 


1014.417 


152.595, 


132.557 




ro77 


-0.430 


925.865 


70.137, 


49.421 




ro77 


-0.398 


1116.000 


36.063, 


36.063 


2.130 


gr88 


-0.383 


833.625 


215.184, 


215.183 


24.798 


be82 


-0.298 


1020.239 


204.589, 


185.769 




ro77 


-0.271 


1117.862 


77.449, 


85.302 




ro77 


-0.222 


1143.000 


61.238, 


61.239 


2.182 


gr88 


-0.144 


1711.456 


381.828, 


319.545 




ro77 


-0.123 


1515.126 


155.710, 


115.090 




ro77 


-0.071 


1278.000 


81.036, 


81.037 


2.440 


gr88 


0.024 


1201.945 


172.364, 


118.204 




ro77 


0.051 


1449.000 


118.832, 


118.833 


2.766 


gr88 


0.138 


1342.800 


147.622, 


147.623 


2.563 


gr88 


0.167 


1515.126 


227.471, 


190.644 




ro77 


0.243 


1234.800 


135.021, 


135.021 


2.357 


gr88 


0.326 


1850.918 


189.560, 


181.435 




ro77 


0.398 


1008.000 


165.611, 


165.612 


1.924 


gr88 


0.474 


1513.772 


227.473, 


213.932 




ro77 


0.602 


1008.000 


165.611, 


165.612 


1.924 


gr88 


0.632 


1171.481 


278.653, 


264.436 




ro77 


0.769 


1291.174 


345.812, 


339.041 




ro77 


0.875 


873.000 


205.206, 


205.207 


1.666 


gr88 


0.923 


1493.462 


471.192, 


494.210 




ro77 


1.081 


969.058 


118.882, 


90.041 




ro77 


1.223 


543.089 


101.008, 85.167 




ro77 


1.377 


647.889 


153.950, 111.570 




ro77 
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Table 3: continued 

log(S[Jy]) norm. cts. + , - sampl. err. ref. 

1.525 768.124 205.943, 169.792 - Toff 

Reference codes. be82:lBenn et al.l dl982t) ; er88: lGruefl Il988f) : 
ro77: iRobertsod |1977l). 



Table 4: Euclidean normalized differential source counts at 610 MHz. 



log(S[Jy]) 


norm. cts. 


+ 


_ 


sampl. err. 


ref. 


-4.252 


4.540 


1.034, 


1.034 


0.805 


ib09 


-4.056 


6.920 


0.780, 


0.780 


0.637 


ib09 


-3.857 


8.790 


0.745, 


0.745 


0.586 


ib09 


-3.662 


11.350 


0.875, 


0.875 


0.646 


ib09 


-3.480 


7.600 


0.559, 


0.559 


0.251 


ga08 


-3.467 


10.330 


0.938, 


0.938 


0.564 


ib09 


-3.432 


11.060 


1.103, 


1.103 


0.593 


bo07 


-3.334 


8.800 


0.578, 


0.578 


0.290 


ga08 


-3.310 


15.100 


2.735, 


2.735 


0.849 


mo07 


-3.271 


10.520 


1.212, 


1.212 


0.569 


ib09 


-3.260 


9.720 


1.137, 


1.137 


0.521 


bo07 


-3.175 


10.100 


0.686, 


0.686 


0.333 


ga08 


-3.081 


12.230 


1.665, 


1.665 


0.656 


bo07 


-3.076 


11.170 


1.710, 


1.710 


0.604 


ib09 


-3.019 


10.300 


0.869, 


0.869 


0.340 


ga08 


-3.013 


15.600 


2.093, 


2.093 


0.877 


mo07 


-2.907 


11.590 


2.113, 


2.113 


0.621 


bo07 


-2.880 


12.720 


2.583, 


2.583 


0.688 


ib09 


-2.865 


12.600 


1.176, 


1.176 


0.415 


ga08 


-2.730 


20.640 


3.814, 


3.814 


1.107 


bo07 


-2.714 


20.200 


3.396, 


3.396 


1.136 


mo07 


-2.710 


14.500 


1.670, 


1.670 


0.478 


ga08 


-2.684 


16.020 


4.161, 


4.161 


0.866 


ib09 


-2.554 


22.520 


5.309, 


5.309 


1.207 


bo07 


-2.554 


21.500 


2.599, 


2.599 


0.709 


ga08 


-2.489 


24.600 


7.371, 


7.371 


1.330 


ib09 


-2.413 


34.400 


7.166, 


7.166 


1.935 


mo07 


-2.391 


27.700 


3.908, 


3.908 


0.913 


ga08 


-2.378 


34.840 


8.908, 


8.908 


1.868 


bo07 


-2.293 


40.230 


13.467, 


13.467 


2.176 


ib09 


-2.243 


28.200 


4.987, 


4.987 


0.930 


ga08 


-2.202 


44.000 


13.478, 


13.478 


2.359 


bo07 


-2.114 


63.600 


15.421, 


15.421 


3.577 


mo07 


-2.098 


57.920 


23.470, 


23.470 


3.132 


ib09 


-2.075 


43.900 


8.425, 


8.425 


1.447 


ga08 


-2.026 


29.390 


14.784, 


14.784 


1.576 


bo07 


-1.937 


78.600 


10.802, 


10.802 


2.537 


ga08 


-1.850 


80.990 


33.354, 


33.354 


4.343 


bo07 


-1.802 


86.100 


14.273, 


14.273 


2.779 


ga08 


-1.678 


168.000 


24.703, 


24.703 


5.423 


ga08 


-1.674 


148.800 


61.272, 


61.272 


7.978 


bo07 


-1.664 


86.500 


26.451, 


26.451 


4.865 


mo07 


-1.552 


180.000 


36.000, 


36.000 


0.000 


ka79 


-1.544 


198.000 


33.613, 


33.613 


6.392 


ga08 


-1.498 


91.120 


64.615, 


64.615 


4.886 


bo07 


-1.418 


202.000 


41.713, 


41.713 


6.521 


ga08 


-1.362 


306.000 


45.000, 


45.000 


0.000 


ka79 
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Table 4: continued 



z(S\Jv]) 


norm. cts. 




+ 




sampl. err. 


ref. 


-1.322 


334.790 


168 


360, 


168.360 


17.951 


bo07 


-1.263 


137.000 


45 


913, 


45.914 


4.423 


ga08 


-1.171 


360.000 


54 


000, 


54.000 


0.000 


ka79 


-1.166 


217.000 


65 


874, 


65.873 


7.005 


ga08 


-1.063 


562.000 


189 


653, 


189.653 


31.609 


mo07 


-1.023 


133.000 


66 


838, 


66.838 


4.293 


ga08 


-0.981 


360.000 


63 


000, 


63.000 


0.000 


ka79 


-0.915 


276.000 


113 


351, 


113.351 


8.910 


ga08 


-0.790 


495.000 


99 


000, 


99.000 


0.000 


ka79 


-0.757 


339.000 


170 


352, 


170.352 


10.943 


ga08 


-0.600 


675.000 


153 


000, 


153.000 


0.000 


ka79 


-0.410 


684.000 


216 


000, 


216.000 


0.000 


ka79 


-0.220 


1314.000 


414 


000, 


414.000 


0.000 


ka79 


0.096 


837.000 


369 


000, 


369.000 


0.000 


ka79 



Reference codes. bo07: 



j2008fl; ib0 9 : llbar et al 
iMoss ct al. (2007]). 



Bondi et all J2007f); ga 8: iGarn et al.l 
' (120091 1; ka79: iKatgertl J1979D ; mo07: 



Table 5: Euclidean normalized differential source counts at 1.4 GHz. 



log(5[Jyj) 


norm. cts. 


+ 




sampl. err. 


ref. 


-4.770 


6.480 


1.510, 


1.510 


0.011 


ow08 


-4.668 


6.500 


0.890, 


0.890 


0.011 


ow08 


-4.602 


3.210 


0.752, 


0.752 


0.676 


ib09 


-4.523 


6.110 


0.580, 


0.580 


0.010 


ow08 


-4.467 


2.600 


1.624, 


1.624 


0.278 


fo06 


-4.398 


4.140 


0.529, 


0.529 


0.471 


ib09 


-4.373 


4.100 


2.048, 


2.048 


0.439 


fo06 


-4.347 


5.750 


0.690, 


0.690 


0.010 


ow08 


-4.340 


0.980 


0.278, 


0.278 


0.171 


se08A 


-4.284 


2.300 


0.688, 


0.688 


0.314 


riOO 


-4.272 


5.400 


2.372, 


2.372 


0.578 


fo06 


-4.244 


2.490 


0.509, 


0.509 


0.136 


ho03 


-4.208 


5.120 


0.511, 


0.511 


0.446 


ib09 


-4.206 


2.580 


0.528, 


0.528 


0.276 


ke08 


-4.180 


3.287 


0.245, 


0.245 


0.176 


bo08 


-4.173 


5.700 


2.476, 


2.476 


0.610 


fo06 


-4.155 


2.900 


0.803, 


0.803 


0.396 


riOO 


-4.151 


1.110 


0.285, 


0.285 


0.193 


se08A 


-4.143 


3.270 


0.447, 


0.447 


0.179 


ho03 


-4.114 


7.010 


0.740, 


0.740 


0.012 


ow08 


-4.092 


4.454 


0.318, 


0.318 


0.239 


bo08 


-4.077 


7.000 


2.706, 


2.706 


0.750 


fo06 


-4.066 


3.660 


0.614, 


0.614 


0.201 


ho03 


-4.060 


3.100 


0.763, 


0.763 


0.423 


riOO 


-4.021 


2.520 


0.466, 


0.466 


0.270 


ke08 


-4.009 


5.060 


0.473, 


0.473 


0.374 


ib09 


-4.000 


4.434 


0.331, 


0.331 


0.238 


bo08 


-4.000 


3.910 


1.255, 


1.255 


0.110 


ho03 


-3.971 


2.100 


0.839, 


0.839 


0.254 


mi85 


-3.960 


5.200 


2.366, 


2.366 


0.557 


fo06 


-3.955 


2.700 


0.727, 


0.727 


0.368 


riOO 


-3.927 


1.610 


0.369, 


0.369 


0.280 


se08A 


-3.914 


4.250 


0.908, 


0.908 


0.120 


ho03 
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Table 5: continued 



gfSfJvl) 


norm. cts. 


+ 




sampl. err. 


ref. 


-3.914 


4.944 


0.371, 


0.371 


0.265 


bo08 


-3.903 


6.240 


1.110, 


1.110 


0.011 


ow08 


-3.836 


4.440 


0.429, 


0.429 


0.125 


ho03 


-3.827 


2.880 


0.546, 


0.546 


0.308 


kc08 


-3.824 


5.147 


0.438, 


0.438 


0.276 


bo08 


-3.810 


4.300 


0.460, 


0.460 


0.298 


ib09 


-3.807 


2.400 


0.698, 


0.698 


0.327 


riOO 


-3.784 


5.100 


2.364, 


2.364 


0.546 


fo06 


-3.783 


3.500 


0.905, 


0.905 


0.423 


mi85 


-3.770 


3.600 


0.404, 


0.404 


0.105 


ci99 


-3.764 


4.960 


0.500, 


0.500 


0.140 


ho03 


-3.738 


5.014 


0.499, 


0.499 


0.269 


bo08 


-3.686 


4.730 


0.412, 


0.412 


0.133 


ho03 


-3.669 


1.820 


0.487, 


0.487 


0.317 


se08A 


-3.650 


4.020 


0.517, 


0.517 


0.216 


bo08 


-3.648 


4.460 


0.940, 


0.940 


0.008 


ow08 


-3.629 


4.400 


1.313, 


1.313 


0.532 


mi85 


-3.627 


2.900 


0.745, 


0.745 


0.396 


riOO 


-3.611 


3.880 


0.549, 


0.549 


0.266 


ib09 


-3.597 


4.670 


0.412, 


0.412 


0.131 


ho03 


-3.569 


5.300 


1.012, 


1.012 


0.158 


er99 


-3.561 


5.300 


0.691, 


0.691 


0.284 


bo08 


-3.543 


4.230 


0.817, 


0.817 


0.453 


kc08 


-3.509 


3.800 


1.285, 


1.285 


0.460 


mi85 


-3.495 


4.610 


0.365, 


0.365 


0.134 


ci99 


-3.487 


4.370 


0.380, 


0.380 


0.123 


ho03 


-3.472 


4.490 


0.712, 


0.712 


0.241 


bo08 


-3.413 


5.480 


0.902, 


0.902 


0.376 


ib09 


-3.384 


4.470 


0.816, 


0.816 


0.240 


bo08 


-3.378 


6.200 


2.587, 


2.587 


0.664 


fo06 


-3.377 


2.530 


0.777, 


0.777 


0.441 


sc08A 


-3.350 


4.640 


0.411, 


0.411 


0.131 


ho03 


-3.347 


3.500 


0.853, 


0.853 


0.477 


riOO 


-3.323 


6.500 


1.519, 


1.519 


0.786 


mi85 


-3.319 


7.700 


0.642, 


0.642 


0.229 


er99 


-3.297 


5.140 


1.008, 


1.008 


0.276 


bo08 


-3.252 


4.780 


0.414, 


0.414 


0.139 


ci99 


-3.215 


4.700 


1.185, 


1.185 


0.322 


ib09 


-3.208 


5.970 


1.319, 


1.319 


0.320 


bo08 


-3.197 


6.150 


0.548, 


0.548 


0.173 


ho03 


-3.120 


4.050 


1.190, 


1.190 


0.217 


bo08 


-3.111 


3.930 


1.275, 


1.081 


0.211 


bo08 


-3.097 


5.400 


1.821, 


1.821 


0.653 


mi85 


-3.070 


6.580 


1.666, 


1.666 


0.705 


ke08 


-3.060 


5.700 


0.528, 


0.528 


0.170 


gr99 


-3.046 


4.230 


1.260 


1.260 


0.007 


ow08 


-3.032 


5.490 


1.607| 


1.607 


0.294 


bo08 


-3.031 


5.483 


1.646, 


1.508 


0.294 


bo08 


-3.025 


7.980 


0.697, 


0.697 


0.225 


ho03 


-3.017 


7.720 


2.175, 


2.175 


0.530 


ib09 


-3.000 


3.550 


0.040, 


0.040 


0.000 


wi97 


-2.991 


6.560 


0.611, 


0.611 


0.191 


ci99 


-2.951 


3.230 


0.050, 


0.050 


0.000 


wi97 


-2.900 


5.190 


0.060, 


0.060 


0.000 


wi97 


-2.889 


11.490 


1.705, 


2.035 


0.616 


bo08 


-2.851 


7.230 


0.080, 


0.080 


0.000 


wi97 


-2.824 


9.600 


3.217, 


3.217 


1.161 


mi85 
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Table 5: continued 





norm. cts. 


-)- 




Seimpl. err. 


ref . 


-2.820 


12.330 


3.922 


3.922 


0.846 


ib09 


-2.804 


7.200 


0.828 


0.828 


0.214 


st99 


-2.801 


7.640 


0.090 


0.090 


0.000 


wi97 


-2.783 


11.500 


0.994 


0.994 


0.324 


ho03 


-2.750 


9.120 


0.110 


0.110 


0.000 


wi97 


-2.738 


11.680 


1.151 


1.151 


0.340 


ci99 


-2.699 


10.040 


0.120 


0.120 


0.000 


wi97 


-2.684 


8.869 


2.457 


2.321 


0.476 


bo08 


-2.650 


12.270 


0.150 


0.150 


0.000 


wi97 


-2.622 


10.630 


5.677, 


5.677 


0.729 


ib09 


-2.600 


12.680 


0.170 


0.170 


0.000 


wi97 


-2.550 


14.160 


0.190 


0.190 


0.000 


wi97 


-2.550 


10.100 


1.432 


1.432 


0.300 




-2.500 


15.410 


0.220 


0.220 


0.000 


wi97 


-2.481 


15.860 


2.033 


2.033 


0.461 


ci99 


-2.479 


14.080 


4.228 


4.061 


0.755 


bo08 


-2.450 


17.100 


0.250 


0.250 


0.000 


wi97 


-2.400 


18.990 


0.290 


0.290 


0.000 


wi97 


-2.350 


20.540 


0.330 


0.330 


0.000 


wi97 


-2.309 


23 980 

— ■ ' I - • .' 1..' V ) 


2.073 


2.073 


0.675 


ho03 


-2.300 


22.890 


0.370 


0.370 


0.000 


wi97 


-2.295 


20.900 


3.260 


3.260 


0.622 


st99 
& 1 yy 


-2.284 


55 000 


14.272 


14.272 


5.890 


kc08 


-2.283 


31.190 


8.643 


8.987 


1.672 


bo08 


-2.250 


25.260 


0.430 


0.430 




wi97 


-2.227 


31.610 


4.436 


4.436 


0.919 


ci99 


-2.200 


27.280 


0.490 


0.490 




wi97 


-2.194 


14.300 


3.820 


3.820 


1.531 


fo06 


-2.150 


31.380 


0.570 


0.570 




wi97 


-2.100 


34.540 


0.650 


0.650 




wi97 


-2.078 


32.370 


12.085 


13.899 


1.736 


bo08 


-2.050 


37.240 


0.740 


0.740 




wi97 


-2.040 


30.000 


6.066 


6.066 


0.892 




-2.000 


42.040 


0.850 


0.850 




wi97 


-1.971 


61.870 


9.600 


9.600 


1.799 


ci99 


-1.950 


46.010 


0.970 


0.970 




wi97 


-1.900 


50.670 


1.110 


1.110 




wi97 


-1.873 


84.690 


31.638 


29.689 


4.541 


bo08 


-1.850 


54.590 


1.260 


1.260 




wi97 


-1.800 


59.780 


1.430 


1.430 




wi97 


-1.783 


52 200 


12.398 


12.398 


1.553 




-1.750 


66.180 


1.650 


1.650 




wi97 


-1.716 


59.070 


14.433 


14.433 


1.718 


ci99 


-1.700 


77.710 


1.940 


1.940 




wi97 


-1.668 


86.270 


41.191 


43.000 


4.626 


bo08 


-1.650 


85.570 


2.220, 


2.220 




wi97 


-1.600 


80.120 


2.350, 


2.350 




wi97 


-1.550 


100.880 


2.870, 


2.870 




wi97 


-1.529 


105.000 


27.279, 


27.279 


3.123 


gr99 


-1.500 


114.910 


3.340, 


3.340 




wi97 


-1.478 


84.390 


11.429, 


11.429 


2.374 


ho03 


-1.461 


109.100 


30.416, 


30.416 


3.172 


ci99 


-1.450 


102.910 


3.450, 


3.450 




wi97 


-1.400 


134.730 


4.300, 


4.300 




wi97 


-1.350 


119.970 


4.420, 


4.420 




wi97 


-1.300 


138.130 


5.170, 


5.170 




wi97 


-1.274 


84.500 


37.883, 


37.883 


2.514 


gr99 
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Table 5: continued 





norm. cts. 


-|- 




sampl. err. 


ref . 


-1.250 


160.170 


6.070 


6.070 




wi97 


-1.206 


162.100 


57.503 


57.504 


4.713 


ci99 


-1.200 


179.240 


7.000 


7.000 




wi97 


-1.150 


172.490 


7.490 


7.490 




wi97 


-1.100 


208 600 


8.980 


8.980 




wi97 


-1.050 


230 730 


10.290 


10.290 




wi97 


-1.019 


122.400 


70.794 


70.794 


3.641 


st99 


-1.000 


223 840 


11.050 


11.050 




wi97 


-0.950 


244.700 


109 631 


109 631 


7.115 


ci99 


-0.950 


231 560 


12.260 


12.260 




wi97 


-0.900 


258 500 


14.120 


14.120 




wi97 


-0.850 


246.720 


15.030 


15.030 




wi97 


-0.800 


328 630 


18.920 


18.920 




wi97 


-0.763 


197.100 


1 39 523 


1 39 523 

-L kj C • *-J *-J 


5.863 




-0.750 


960 160 


18.350 


18.350 




wi97 


-0.700 


308 060 


21.770 


21.770 




wi97 


-0.695 


236 400 


167.241 


167.241 


6.874 


ci99 


-0.650 


321 100 


24.230 


24.230 




wi97 


-0.600 


371 500 


28.410 


28.410 




wi97 


-0.550 


375 950 


31.160 


31.160 




wi97 


-0.500 


305 91 


30.640 


30.640 




wi97 


-0.450 


296 630 


32.890 


32.890 




wi97 


-0.400 


340 950 


38.440 


38.440 




wi97 


-0.350 


351 930 

Jul . J Jv 


42.580 


42.580 




wi97 


-0.300 


364.210 


47.220 


47.220 




wi97 


-0.250 


328 270 


48.880 


48.880 




wi97 


-0.200 


632 160 


73.940 


73.940 




wi97 


-0.150 


482 540 


70.430 


70.430 




wi97 


-0.100 


273 390 


57.790 


57.790 




wi97 


-0.050 


490 070 


84.350 


84.350 




wi97 


0.002 


286 500 


65 800 


80 300 

L I KJ . t 1 \ ) \ J 




wi97 


0.050 


339 200 


81.900 


92.700 




wi97 


0.100 


286 400 


75.900 


97.000 




wi97 


0.149 


512 800 


111.800 


1 36 700 




wi97 


0.198 


317 500 


115.400 


1 07 500 




wi97 


0.246 


291 800 

-./ 1 .LI L' V / 


1 33 000 


107 600 




wi97 


0.296 


436 900 

it' L' . - / L' V 7 


181.700 


1 39 400 




wi97 


0.350 


31 8 600 


29.375 


31.090 


9.191 


br72 


0.351 


283.600 


114.100, 


161.800 




wi97 


0.449 


851.100 


276.900, 


287.400 




wi97 


0.452 


254.925 


23.457, 


18.821 


7.354 


br72 


0.498 


568.200 


213.800, 


289.700 




wi97 


0.547 


488.800 


300.500, 


246.900 




wi97 


0.548 


212.512 


32.643, 


30.965 


6.130 


br72 


0.596 


311.300 


295.300, 


214.080 




wi97 


0.700 


332.400 


268.400, 


299.410 




wi97 


0.704 


206.865 


31.751, 


32.569 


5.967 


br72 


0.798 


557.200 


646.800, 


362.300 




wi97 


0.891 


141.885 


23.999, 


22.338 


4.093 


br72 


1.359 


150.570 


27.903, 


25.461 


4.343 


br72 



Refere nce c odes. bo08: iBondi et alj ll2008f): br72:lBridle et alJ 
119721); ci99:ICilieei et al.l dl999t ); fo06: lFomalont et al l d2006f) ; 
gr99: iGruppioni et all lll999al); ho03: iHopkins et al.M2003l) ; 
ib09: llbar et al.l j2q09T); ke 08: iKellermann et al.l (120081); mi85 : 
iMitchell fc Condon! j 1985ft ; ow 08: lOwen fc MorrisoJ <2008ft ; 
riOO: iRichardsl [120001 1; se08A: ISevmour et al.l l|2008n : wi97: 
IWhite et al.l <1997f>. 
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Table 6: Euclidean normalized differential source counts at 4.8 GHz. 





norm. cts. 


+ 




sampl. err. 


ref. 


-4.750 


0.810 


0.562, 


0.562 


0.156 


fo91 


-4.587 


1.080 


0.341, 


0.341 


0.208 


fo91 


-4.347 


1.890 


0.453, 


0.453 


0.364 


fo91 


-3.928 


1.161 


0.455, 


0.455 


0.149 


do87 


-3.870 


0.990 


0.331, 


0.331 


0.191 


fo91 


-3.862 


1.975 


1.192, 


0.977 


0.381 


fo84 


-3.682 


1.503 


0.416, 


0.416 


0.193 


do87 


-3.483 


2.136 


0.967, 


0.971 


0.412 


fo84 


-3.319 


3.807 


1.116, 


1.116 


0.489 


do87 


-2.919 


4.383 


2.161, 


2.134 


0.845 


fo84 


-2.886 


9.027 


2.823, 


2.823 


1.158 


do87 


-2.857 


6.750 


3.160, 


3.160 


1.301 


fo91 


-2.721 


6.210 


1.726, 


1.726 


0.235 


wr90 


-2.538 


8.370 


2.987, 


2.987 


0.317 


wr90 


-2.352 


9.630 


3.618, 


3.618 


0.364 


wr90 


-2.137 


14.850 


6.235, 


6.235 


0.562 


wr90 


-1.961 


29.700 


11.037, 


11.037 


1.124 


wr90 


-1.921 


20.340 


6.982, 


6.982 


0.285 


pa80 


-1.796 




8.675, 


8.675 


0.394 


pa80 


-1.785 


34.151 


2.586, 


2.586 


0.243 


al86 


-1.668 


40.000 


2.003, 


2.003 


0.104 


gr96 


-1.668 


44.165 


4.130, 


4.130 


0.314 


al86 


-1.648 


46.350 


7.216, 


7.216 


0.650 


pa80 


-1.536 


40.520 


4.083, 


4.083 


0.288 


al86 


-1.523 


45.000 


2.001, 


2.001 


0.046 


gr96 


-1.389 


60.322 


6.204, 


6.204 


0.429 


al86 


-1.387 


54.270 


8.814, 


8.814 


0.761 


pa80 


-1.372 


55.000 


3.001, 


3.001 


0.057 


gr96 


-1.222 


64.000 


3.000, 


3.001 


0.056 


gr96 


-1.206 


63.551 


7.664, 


7.664 


0.452 


al86 


-1.071 


70.000 


4.000, 


4.000 


0.061 


gr96 


-1.020 


73.650 


12.119, 


12.119 


0.524 


al86 


-0.903 


'vo r\r\r\ 

78.000 


4.001, 


4.001 


0.068 


gr96 


-0.894 


75.960 


15.542, 


15.542 


1.066 


pa80 


-0.757 


81.000 


5.000, 


5.001 


0.070 


gr96 


-0.602 


84.000 


5.001, 


5.001 


0.073 


gr96 


-0.398 


90.000 


5.001, 


5.001 


0.078 


gr96 


-0.222 


84.000 


6.000, 


6.000 


0.073 


gr96 


-0.071 


107.000 


9.001, 


9.000 


0.093 


gr96 


0.032 


103.410 


10.710 


8.100 




ku81 


0.097 


96.000 


9.000, 


9.000 


0.083 


gr96 


0.102 


92.250 


6.210, 


7.263 




ku81 


0.173 


84.987 


5.733, 


5.400 




ku81 


0.243 


83.000 


12.000, 


12.000 


0.072 


gr96 


0.278 


79.587 


9.675, 


8.622 




ku81 


0.398 


62.000 


11.000, 


11.000 


0.054 


gr96 


0.412 


87.813 


9.117, 


12.042 




ku81 


0.539 


79.587 


11.133, 


16.317 




ku81 


0.764 


62.244 


11.088, 


12.762 




ku81 



Refere nce cod es. al86:lAltschulerl dl986fl : do87: Donnel ly et alJ 
J1987T) ; fo84 : iFomalont et al.l lll984ft ; f o91: IFomalont et al.l 
jl99ll i : gr96:lGregorv et al.l 1 1996 ): ku81: lKuehr et al.l (Il98lh: 
pa80: iPaulinv-Toth et al.l l|1980li : wr90: IWrobel fc Krause) 
dl990D. 
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Table 7: Euclidean normalized differential source counts at 8.4 GHz. 



1 or T, T 1 "\ 


norm. cts. 


+ 




sampl. err. 


ref. 


-4.863 


0.231 


0.050, 


0.050 




fo02 


-4.717 


0.520 


0.385, 


0.385 


0.106 


wi93 


-4.611 


0.192 


0.043, 


0.043 




fo02 


-4.480 


0.530 


0.263, 


0.263 


0.108 


wi93 


-4.348 


0.288 


0.064, 


0.064 




fo02 


-4.192 


1.320 


0.577, 


0.577 


0.270 


wi93 


-3.710 


0.441 


0.096, 


0.096 




fo02 


-3.678 


0.834 


0.352, 


0.352 


0.089 


he05 


-3.481 


1.190 


0.547, 


0.547 


0.244 


wi93 


-2.745 


2.787 


1.125, 


1.125 


0.296 


he05 


-2.699 


8.716 


3.107, 


3.107 


0.393 


he05 


-2.474 


20.374 


11.961, 


11.961 


2.165 


he05 


-2.360 


10. 123 


4.550, 


4.550 


0.457 


he05 


-2.160 


10.809 


5.426, 


5.426 


0.488 


he05 


-2.034 


20.953 


9.918, 


9.918 


2.227 


he05 


-1.901 


18.039 


8.108, 


8.108 


0.814 


he05 


-1.654 


19.820 


3.769, 


4.360 




wi93 


-1.618 


41.222 


20.695, 


20.695 


1.860 


he05 


-1.457 


40.823 


8.986, 


6.514 




wi93 


-1.348 


92.022 


65.202, 


65.202 


4.153 


he05 


-1.250 


55.132 


15.549, 


13.185 




wi93 


-1.052 


31.982 


6.081, 


5.758 




wi93 


-0.855 


43.182 


8.211, 


9.492 




wi93 


-0.658 


62.820 


11.945, 


11.321 




wi93 


-0.461 


62.951 


19.796, 


17.368 




wi93 


-0.263 


82.912 


23.387, 


19.831 




wi93 


-0.056 


61.674 


9.908, 


2.992 




wi93 


0.122 


63.343 


6.625, 


7.380 




wi93 


0.320 


58.911 


6.161, 


5.566 




wi93 


0.527 


43.823 


8.309, 


7.898 




wi93 


0.726 


24.797 


10.319 


9.330 




wi93 



Reference codes. fo02: iFomalont et al 1 (|2002|) 
Hcnkcl fc Partridge (2005): wi93: Windhorst ct al. 



he05: 



(1993). 



Table 10: Euclidean n ormalized differential source counts from WMAP maps 
l|Massardi et al.ll2009h . 



log(5[Jyj) 


norm. cts. 


+ , - 


v (GHz) 


0.111 


48.993 


3.720, 3.721 


23 


0.333 


37.111 


4.639, 4.639 


23 


0.556 


32.480 


6.371, 6.370 


23 


0.778 


32.297 


9.324, 9.323 


23 


1.000 


11.597 


7.538, 10.594 


23 


0.111 


75.708 


8.316, 8.317 


33 


0.333 


38.878 


4.646, 4.647 


33 


0.556 


28.717 


5.863, 5.862 


33 


0.778 


28.357 


8.550, 8.550 


33 


1.000 


16.662 


9.097, 12.730 


33 


1.222 


23.931 


15.556, 21.861 


33 



continued on the next page 
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Table 10: continued 





norm. cts. 


+ , - 


v (GHz) 


0.111 


58.115 


5.670, 5.670 


41 


0.333 


34.633 


4.364^ 4.363 


41 


0.556 


28.425 


5.802, 5.802 


41 


0.778 


20.413 


7.068, 9.138 


41 


1.000 


21.989 


10.544, 14.474 


41 


1.222 


23.688 


15.397, 21.639 


41 


1.444 


25.516 


21.459, 25.516 


41 


0.333 


26.078 


3.805, 3.804 


61 


0.556 


21.307 


5.022, 5.022 


61 


0.778 


10.201 


4.892, 6.715 


61 


1.000 


16.483 


9.000, 12.593 


61 


1.222 


23.674 


15.388, 21.627 


61 


1.444 


25.502 


21.444, 25.502 


61 
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